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PREFACE 
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Analysis Branch, Earth Resources Research Division, Space and Life Sciences 

Directorate of the Lyndon B. Johnson Space Center, National Aeronautics and 

Space Administration. All work was performed under Contract NAS 9-15800. 

The following scientists and other personnel contributed to this work: 

1. 6. Badhwar of the Lyndon B. Johnson Space Center developed some of the 
techniques employed in the labeling procedure that is provided in this 
document. He is responsible for the multitemporal adaptation of the edge 
gradient technique used in eliminating impure pixels, for the greenness 
over brightness profile technique used in determining crops or noncrops, 
and for the direction of the computer adaptation of the growth stage 
prediction models used in predicting growth stages. Successful operation 
of the labeling process is highly dependent upon these techniques. 

2. J. G. Carnes of Lockheed Engineering and Management Services Company, 
Inc., provided many constructive ideas for the preparation of this 
document. 

3. G. L. Clouette, R. L. Fagan, C. L. Horton, J. L. Lastres, J. H. McCormac, 
C. A. Sivillo, and M. A. Tompkins of Lockheed Engineering and Management 
Services Company, Inc., provided invaluable computer programming and data 
processing support during the development and testing of this procedure. 

4. B. A. Tolbert of Lockheed Engineering and Management Services Company, 
Inc., participated in testing the integrated labeling and classification 
system and provided many ideas for improving clarity in the procedures 
documentation. 
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1. INTRODUCTION 


1.1 OBJECTIVES 

The purposes of this report are (1) to document some essential processes for 
the development of a green-number-based logic for Identifying (labeling) crops 
In Landsat Imagery and (2) to record the supporting data and subsequent 
conclusions that resulted from development of a specific labeling logic [l.e., 
for corn and soybean crops In the United States (ref. 1)]. 

The corn and soybean labeling logic developed, while possibly being useful to 
other crop classification systems, was designed specifically to Isolate and 
Identify pure pixels to train a specific classifier. The classifier to be 
trained uses profile parameters of temporal greenness curves as the basis for 
separating crops of Interest In Landsat data. In addlvlon to high purity. It 
was desired that the labeling accuracy of crops also be high In order to keep 
the classification errors due to mislabeling at a low level. Earlier experi- 
mentation (ref. 2) had shown th*c this classifier could accurately clissify 
both pure and mixed pixels In •' Landsat scene with good accuracy when trained 
only with "superpure" ground truth pixels. It was now desired to see whether 
comparable accuracy could be attained using training pixels that were both 
Isolated and labeled using either an automated or semlautomated labeler. 

1.2 MAJOR CONSIDERATIONS 

The identification of crops In Landsat imagery by Lockheed personnel 
supporting the Johnson Space Center (JSC) Earth Observations Division [l.e., 
during the Large Area Crop Inventory Experiment (LACIE), the Transition Year 
project, and the Agriculture and Resources Inventory Surveys Through Aerospace 
Remote Sensing (AgRISTARS) program] has always been accomplished by the 
recognition o f 2 temporal pattern of characteristics through the growing 
season of a particular crop of interest. If the time periods of occurrence of 
the recognizable characteristics can be Isolated accurately year to year, and 
these periods do not overlap prohlbltl vely with periods in which similar 
characteristics exist In other crops, the crop of Interest can be consistently 
and accurately Identified. 
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Since 1978, efforts to improve crop labeling have concentrated on methods for 
Increasing objectivity and repeatability, as well as methods for Improving 
utilization of computers. The Kauth transformations of greenness and bright- 
ness (ref. 3) have :ontr1buted much to these efforts. The temporal sequences 
of colors (l.e., noired-red-nonred) signifying nonvegetatlon-veoetatlon- 
nonvegetatlon present In the false-color film products that are used In 
AgRISTARS are sufficiently relatable to plant greenness (as green numbers) and 
brightness to allow Identification of the same sequences using greenness ar 
brightness. Since greenness and brightness are numerical quantities generated 
from the Landsat channel values, efforts have been directed at creating 
computer logic using these numbers to Identify the sequences In an objective 
and repeatable way. Reference 4 documents a labeling procedure for spring 
small grains that employs greenness and brightness In this context, and the 
research addressed In section 3 of this report provides a greenness 
methodology for Identifying corn and soybeans as summer crops using the same 
basic technique. Brightness, in combination with greenness, Is shown In 
section 3.4 to oe valuable in separating corn and soybeans, both from each 
other and from certain other summer crops. 

Of utmost Importance, and perhaps the most elusive to predict, Is the time 
frame of occurrence during l new crop year of the identifying characteristics 
of a crop. While the spectral characteristics of a crop, as recorded by the 
Landsat, seem to have a direct relationship to the phenology of the crop, 
accurate modeling of all the necessary phenologlcal factors over large areas 
is very difficult. Section 3.7 covers the research carried out to adapt two 
growth stage prediction models for the purpose of predicting the key time 
periods necessary to the process of Identifying corn and soybeans. It was 
desired by the National Aeronautics and Space Administration (NASA) technical 
monitor that these models be evaluated for use. The sensitivity of the 
procedure to the accuracy of the model and the compensating corrections are 
discussed. 

The size of the area over which a crop must be Identified is an Important 
consideration both from the standpoint of predicting time of occurrence of 
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Identifiable stages of che crop arid from the standpoint of the number of types 
of confusion crops that must be dealt with. The primary objective of this 
effort was to develop a labeling procedure for the designated U.S. Corn Belt 
States of Iowa, Illinois, and Indiana. The segments In the developmental data 
set were therefore chosen primarily to address the conditions In this area. 

But In order to Increase flexibility, some peripheral areas and one area In 
the Deep Soutn that grows primarily soybeans were considered. Logic was then 
developed for labeling these areas using the selected segments. Other areas 
would have been addressed If data had been available. It Is felt that If 
growth stages can be predicted adequately and the necessary Landsat acquisi- 
tions are available segments would be processable In most areas of the 
United States. 

Finally, the output from a labeling procedure designed to provide training 
Information for a classification algorithm must take Into account the require- 
ments of the classifier. The Ho-Kashyap-type classifier used In profile 
parameter classification requires an equal number of pure dots In each crop 
category to be classified. The technical monitor established 90 percent or 
greater as the labeling accuracy goal for these dots. The steps taken during 
this research effort to Isolate only pure dots are covered In sections 4.2.3 
and 4.2.4 of this document. Numerous factors were considered In arriving at 
methodology to achieve the required labeling accuracy. Selection of -he 
appropriate subset from the available pure dots for use by the classifier was 
accomplished with a random number generator. The number of dots required per 
crop category was establisned by separate research (ref. 2). 
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2. DATA SET 


2.1 BASIC GUIDELINES 

Guidelines from the NASA monitor were that a labeling procedure should be 
developed primarily for labeling 5- by 6-nautlcal-mlle segments In the Corn 
Belt States of Iowa, Illinois, and Indiana; but. If possible, the procedure 
should also be usable In other areas of the United States that grow corn and 
f soybean crops. Data availability for the full growing seasc.i would be 

assumed. Additionally, the landsat data for the years 1 78 and 1979 were 
designated for use In developing the labeling procedure. At least six 
\ segments for the year 1979 were also to be set aside for use In a verification 

test of the end-to-end labeling and proportion estimation procedure that would 
follow the developmental effort. These latter six segments were not to be 
used in developing the labeling procedure. 

Additionally, ground truth data for the Iowa 1980 season were not to be used, 
because Iowa segments were reserved for use in a large-area demonstration of 
the procedure. 

Because of the large-scale demonstration requirement, the procedure would 
therefore have to be executable using 1980 data. 

2.2 DATA SET CHOSEN FOR DEVELOPMENT 

Figure 2-1 is a map of states containing segments used in labeling procedure 
development. Counties In which the segments are located are shaded. Table 2-1 
t Is a tabulation of the segments in Increasing numerical order and shows the 

acquisitions used for generating spectral aids for use in the developmental 
techniques. This table also shows the percentage (oy year) of area devoted to 
corn and soybean acreage that was Inventoried by the ground surveyor. 

Selection criteria for segments Included the following factors: 

1. Geographic distribution throughout Iowa, Illinois, and Indiana for the 
years 1978 and 1979. 
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TABLE 2-1. 'DEVELOPMENTAL SEGMENT LIST 
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l Figure. 2-1.- 


Counties containing developmental segments. 
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2. Cloud-free acquisitions for each segment to bracket the periods of plant- 
ing, peak green-up, and harvest/postharvest for both com and soybeans. 

3. At least one cloud-free acquisition for each segment showing maximum 
interpreted color separation of color and soybeans (when both were 
present) in the production film converter (PFC) product 1. 

4. At least one segment having com but no soybeans. 

5. At least one segment having soybeans but no corn. 

6. The mix of confusion crops would be as typical as possible of the area in 
which the segment was located. 

In addition, it was decided that segments would be included in the peripheral 
areas of Minnesota and Nebraska. The separability of spring small grains, 
sunflowers, and oth>r possible confusion such as hay could thus be examined in 
the U.S. northern Great Plains. 

In order to find a segment that grew soybeans but not corn, it was necessary 
to use Mississippi segments. The additional benefit of examining the 
separability of cotton accrued from using segments in that state. Segment 
data from both years were used in those cases where sufficient acquisitions 
were available. 

2.3 DATA SET CHOSEN FOR VERIFICATION TEST 

At the time the developmental segments were selected, nine segments in Iowa 
and Indiana were reserved for potential use in the verification test, in the 
hope that at least six of them would have sufficient acquisitions within the 
correct time periods for processing with the procedure. However, the 
procedure that resulted from the developmental research had selection criteria 
(see section 3.5) that eliminated all but four of the segments. Enough 1978 
segments were selected at random from a listing, without regarding location, 
to make up the deficit. Table 2-2 lists all the segments, and figure 2-2 is a 
plot of the segments. 


2-4 


TABLE 2-2.- TEST SEGMENT LIST 
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Figure 2-2.- Counties containing test segments. 
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2.4 DATA PREPARA. ION 


In order to establish a basis for relating the Kauth transformation values of 
the corn and soybean crops to Identifiable characteristics of the c^ops In the 
PFC product 1 (Imagery), It was necessary (1) to Identify pure corn and 
soybean dots for use in the greenness and brightness study and (2) to generate 
spectral aids of the dots for studying the greenness characteristics of these 
crops versus confusion (pasture, other crops, and nonagriculture). Study was 

, necessary temporally, using multiple acquisitions and single (but critical) 

acquisition dates. 

i Green number/brightness listings of pixels (dots) at the 209 grid intersec- 

tions of the PFC products and also plots of green numbers through time and 
brightness through time of the dots were first generated. (See appendix A.) 
All 209 dots were then manually labeled with ground truth using ground truth 
annotated photographs to determine the Identity of the fields or areas 
containing the dots. Following this, each dot was evaluated for purity using 
PFC products 1 and 3 and the green number through time trajectories. A pure 
dot was taken to be one that was located entirely within the same field (or 
area type) on all of the selected acquisitions. After the pure dot ground 
truth was determined, labeled scatter plots of corn and soybeans (C and Y) and 
other (N) categories were generated for use In single-date evaluation 
(appendix B) using only these pure dots. 

Labeling procedures presently In use at JSC employ green numbers rather than 
the Kauth greenness directly. Green numbers more closely approximate the 
greenness of the plants themselves and are generally more stable as a gauging 
device from acquisition to acquisition and from segment to segment. Green 
numbers represent total greenness minus bare soil greenness. Occasionally, 
there are errors in the calculation of green numbers, or nonrepresentative 
green numbers result from correct calculations because the formula Is based on 
the premise that there will be at least some bare soil In the segment. In the 
absence of bare soil, some vegetation (of low greenness value) will be 
assigned bare soil greenness, and all of the calculations of green numbers 
will be biased to give values too low. It It therefore necessary to verify 
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the validity of green numbers In each acquisition before basing decisions on 
green numbers. Before using the spectral aids described In the above 
paragraph to derive the labeling decision rules In section 3, the following 
criteria were applied to green numbers for each segment In the data set: 

1. Green numbers of water or clouds should be zero or negative values. 

2. Green numbers of plowed ground should be near zero or slightly negative 
quantities (within 0 ± 5). 

3. Green numbers of greened-up crops and other solid red vegetation should be 
greater than or equal to 13. 

4. Green numbers on the time plots of pure crop dots should form continuous 
normal growth curves [fig. 2-3(a) and 2-3(b)]. If summer crop dots 
throughout the segment show a deviation similar to the pixel In 

figure 2-3 ( b) , a correction In the soil line Is necessary. 

5. The lower summer crop cluster should not appear to be chopped off at the 
bottom in the labeled scatter plot. 

Soil line corrections were estimated for those acquisitions appearing to have 
bad green numbers, and spectral aids were then regenerated prior to using the 
data. 
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Figure 23.- Growth curves for crops. 
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3. RESEARCH 


3.1 GENERAL 

Under the beslc guidelines provided, the responsibility of Lockheed personnel 
was to arrive at an automated or semi automated labeling logic for labeling a 
small number of pure (field) pixels In a Landsat scene with high accuracy. 
These labeled pixels would be used to train a Ho-Kashyap-type classifier that 
would classify pixels based on greenness profile parameters. Based on the 
major considerations discussed In the Introduction and the experience of the 
author, the research would need to address the following areas In order to 
arrive at the basic elements of a procedure that would meet these 
requirements: 

1. A repeatable logic fcr Identifying corn and soybeanr. as summer crops, 
either directly or by first establishing that the fields contain crops. 

2. A repeatable logic for dU '•“tely Identifying corn and/or soybeans, either 
directly or by going throuv,' Lhe step(s) described in subparagraph 1. 

3. The Identification of specific Acquisitions required for identifying the 
crops and a method for selecting them (analyst, normal crop calendar, 
models, or other method). 

4. Screening of acquisitions to eliminate possible errors caused by clouds, 
excessive misregistration, or bad data (banding, skip, etc.). 

The following subsections address the above essential topics. Although 
acquisition screening and selection steps would precede crop logic execution 
In a procedure, the derivation of crop Identification logic Is documented 
first because acquisition requirements are necessarily based on this logic. 

3.2 SELF-IHPOSED ACCURACY GUIDELINES 

Following labeling logic development, It would be necessary to Integrate this 
logic first with any accepted machine programs for pure pixel Isolation; and, 
second, with the existing programs for generating profile parameters (e.g., 
CLASFYG; ref. 5) In order to execute a complete labeling and classification 
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system. However, because of the constraints of the resources allocated and 
the short period scheduled for research and Integration, there would be little 
time available for making trial runs of the completed classification system 
prior to formal testing. It was necessary, therefore, to establish some basic 
accuracy criteria for each labeling step to be researched In order to 
reasonably assure a high degree of labeling success during Initial shakedown 
of the Integrated system: 

1. Acquisitions with more than 5 percent cloud cover, cloud shadows, or other 
bad data would not be used. 

2. Segments with more than 5 percent confusion crops with unpredictable 
spectral characteristics should not be used (l.e., no more than 5 percent 
of the corn and soybean labels could In reality be other crops of 
nonseparabil 1 ty ) . 

3. No more than 10 percent of a category of Interest could be mislabeled 
because of acquisition selection. This takes Into account the effects of 
executing all of the labeling logic using all labeling acquisitions 
selected. 

It was hoped that cumulative errors from all of the above causes would be 
under 10 percent In the majority of the developmental segments. If this 
occurred, the overall guidelines for labeling accuracy established by NASA 
could be met. 

3.3 LOGIC FOR IDENTIFYING CORN AND SOYBEANS AS SUMMER CROPS 

The procedure used during the Transition Year project first separated cropland 
from noncropland by using a color sequence logic executed In a manual mode 
(ref. 6). While this was considered as a viable possibility for use In the 
procedure being developed, It was desired to see whether a green number logic 
could be developed that would Identify corn and soybeans directly as summer 
crops without first identifying the fields as cropland. This logic, to be 
successful, would have to exclude noncropland from being called summer crops. 
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3.3.1 


GREEN NUMBER LOGIC FOR IDENTIFYING CORN 


First, a sequence of green numbers (GN's) for Identifying corn was somewhat 
arbitrarily established In each of three time periods. The author's 
experience with spring small grain Identification and some study of the 
Imagery resulted In the following: 

1. GN < 10 during the planting period; this period would extend from the time 
natural vegetation reached GN > 10 until the time when no more than 

5 percent of corn had reached GN > 10. 

2. GN > 13 during the period of maximum greenness; this period could begin 

when 95 percent of corn reached a green number of 13 and Ideally would be 

past the time when 95 percent of spring small grains have peaked and 

dropped to a green number less than 13; and It could extend to the point 

at which no more than 5 percent of com (while maturing) had dropped In 
greenness to GN < 13. Corn with green numbers less than 13 can give the 
spectral appearance of being trees. 

3. GN < 10 during harvest/postharvest; at least 95 percent of corn must have 
dropped to GN < 10 at the end of the season. 

3.3.2 GREEN NUMBER LOGIC FOR IDENTIFYING SOYBEANS 

Second, a sequence of green numbers for Identifying soybeans was arbitrarily 
established: 

1. GN < 10 during the planting period; this period would extend from the time 

natural vegetation had reached GN > 10 until the time when no more than 

5 percent of soybeans had reached GN > 10. 

2. GN > 13 during the period of maximum greenness; this period could begin at 

a point where at least 95 percent of the soybeans had reached GN > 13, and 

It could extend to the point at which no more than 5 percent of soybans 
had again dropped to GN < 13. 

3. GN < 10 during harvest/postharvest; at least 95 percent of soybeans would 
have dropped to GN < 10 at the end of the season. 
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The green number thresholds established were subject to adjustment at a later 
point In proceu- r e development. If necessary. Some basis was needed at this 
point for examining the data at critical time periods, to obtain an estimate 
of the length of the critical periods and to make some correlations with other 
data, as explained In section 3.3.4. 

3.3.3 DETERMINATION OF LENGTH OF CRITICAL PERIODS 

The maximum duration o f the critical periods for the two crops was determined 
using the following methods: 

1. Planting periods: With a random selection of ground truth dots for each 

crop, the percentage accuracy of crop Identification was determined using 
each acquisition that bracketed the period and during the period (using a 
green number threshold of 10). Particular attention was paid to 
acquisitions giving marginal accuracy 'In order to determine approximate 
limits of the period. Evaluation of the data for both 1978 and 1979 
Indicated that a period of 45 days for corn and a period of 30 days for 
soybeans met the specified conditions. 

2 . Period of maximum greenness: For each of the randomly selected dots (15 

where possible), the greenness through time plot for the segment was used 
to approximate the date of maximum greenness. The date of occurrence of 
maximum greenness for each crop In eacN segment was taken to be the mean 
of the dates determined for the 15 dots. The accuracy of Identification 
of corn (using the green number thr-shold of 13) was determined for each 
acquisition that bracketed *he period of peak greenness and during the 
period of Deak greenness. The same accuracies were determined for 
soybeans using the appropriate acquisitions. Particular note was taken of 
acquisitions with marginal accuracies In ordt-r to approximate the limits 
of this critical period for each crop. Evaluation of the dots for both 
1978 and 1979 Indicated that a period of about 30 days for corn and a 
period of about 25 days for soybeans met the specified conditions. 

3. Harvest/postharvest period: For each acquisition after crop maturity, and 

using the 15 random dots of each crop, the accuracy of crop identification 
was determined using the green number threshold of 10. The length of this 
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period Is of Indefinite length after harvest, but a length of 35 days was 
arbitrarily chosen to Illustrate the advantages and disadvantages of using 
this time period In a temporal Identification logic. 

3.3.4 GRAPHIC DISPLAY AND EVALUATION OF THE DATA 

As part of the evaluation of whether the defined periods could be used In an 
automated or semlautomated procedure for Identifying corn and soybeans as 
simmer crops, It was necessary (1) to determine In which order to apply the 
green number criteria of the periods In order to minimize omission and 
commission errors, (2) to determine the relationships of the time periods for 
the two crops to each other and decide whether there were conditions under 
which the periods for one crop could be used to Identify both crops, and 
(3) tc determine how the periods could be related to standard growth stages on 
a year-to-year basis In order to help predict and determine the time of 
occurrence of the periods during each new crop year. 

3. 3.4. I The Order In Which to Apply the Green Number Thresholds 
For FhV Tlme Periods 

It was desired to keep the number of required acquisitions for labeling as low 
as possible and still meet the labeling accuracy requirement. This serves to 
minimize the number of segments that will be declared unprocessable, because 
certain required labeling acquisitions are missing. 

App',ir>j the green number thresholds In two periods and In the following order 
gen « -My will result In the least omission of corn and soybeans as summer 
crops wh n using one acquisition per period: 

Harvest Maximum green 

However, becausp r t confusion crops and some trees and pasture will have a 
green number of 10 or less during the harvest period, a large number of these 
categories will be called summer crops (committed) and will not be eliminated 
by the acquisition In the maximum green period for corn and soybeans. One or 
more preemergence acquisitions will be required t.u eliminate the nonsummer 
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dots thus committed, requiring a total of three labeling acquisitions as a 
minimum. Two examples of this are segments 0141 and 0809 for 1978: 


Harv. acq. 
Segment commission 


After max. 
green acq. 


After preemerg. 

acq. Omission 


0141 31 dots 22 dots 2 dots 

0809 27 dots 17 dots — 


1 dot 
1 dot 


In the case of segment 0141, the available preemergence acquisition largely 
corrected for the high commission rate, whereas in segment 0809 none was 
available. Use of a harvest acquisition (although minimizing omission) has 
the additional disadvantage of not allowing the labeling and classification 
system to be utilized until the crops have been harvested. 

Use of the following combination with a single acquisition in each period 
results in the highest omission, but results in low commission of nonsummer 
crops: 

Maximum green Preemergence 

The maximum greenness acquisition picks up almost all of the corn and soybeans 
as summer crops and also eliminates a high percentage of nonsummer crops as 
"N". The preemergence acquisition then eliminates virtually all of the 
remaining nonsummer dots, but also calls some corn and soybean dots nonsunmer 
crops. The omission problem was found to be particularly true in states other 
than Nebraska, Iowa, and Illinois. Using two or three preemergence acquisi- 
tions and applying the logic "is the green number 10 or less than 10 on any 
acquisition" was found to reduce omission to acceptable levels, because this 
increases the chance of detecting the planting signature of the crops. 

Segment 0843 for 1978 (in Indiana) is an example of the two-acquisition case: 



Omission 

Omission 

Omission 


after max. 

after one 

combining two 

Segment 

green acq. 

preemerg. acq. 

preemerg. acqs. 

0843 

1 dot 

7 dots 

2 dots 
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It was decided that the best order for combining acquisitions in the logic was 
maximum greenness acquisitions followed by preemergence acquisitions, because 
a total of two acquisitions in Nebraska, Iowa, and Illinois would be 
sufficient to minimize omission and commission; and three or more acquisitions 
would only be required in areas outside these states. 

3. 3. 4. 2 Relationships of the Time Periods for Corn and Soybeans to Each 
Other; Common Application 

In order to display the data used to define the critical time period lengths, 
to show the relationships of the times of occurrence (overlap, gap, etc.) of 
the periods for the two crops, and to establish a basis for relating the 
periods to standard qrowth stages, the graphics in appendix C were prepared. 
These graphics are for corn and soybeans for 1978 and 1979. Each basic 
graphic has recorded on it: 

1. The segment number and location (county, state). 

2. The crop growth stages nearest the time of occurrence of the critical 
periods defined earlier. Dates of these stages are to be taken as points 
of reference for plotting the segment acquisitions on a relative time 
scale. The reference stages are from the best available CRD or state 
nominal crop calendar for the applicable y e, r that most nearly applied to 
each segment. Reference stages for the periods were: 


Crop 

Period 

Reference stage(s) 

Corn 

Preemeroenca 

50% planted 

Corn 

Maximum greenness 

50% tasseled, 50% dented 

Corn 

Harvest 

50% mature 

Soybeans 

Preemergence 

50% planted 

Soybeans 

Maximum greenness 

50% podded, 50% turned 

Soybeans 

Harvest 

5U% mature 


3. Each acquisition date for each segment plotted in relation to the refer- 
ence growth stages. 
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4. The date of occurrence of growth stages recorded by the Agricultural 
Stabilization and Conservation Service (ASCS) enumerator during the ground 
data surveys. These dates were not always available. 

5. Comments addressing the adequacy of acquisitions at early and late 
extremes of the critical periods. 

6. Approximations of the limits of the critical periods. 

7. References to the crop calendars used. 

Information was added to the graphics incrementally as research progressed. 
Reference is made to them in later sections of this paper when they apply. 

Inspection of the reference nominal crop calendars and the green number 
behavior of the two crops shows that in all areas where segments were selected 
soybeans were planted later than corn and peaked in greenness later than corn. 
Corn, however, remained in the field later than soybeans. There was signifi- 
cant overlap of the critical periods in the U.S. Com Belt states and 
Minnesota. The approximate limits of the critical corn periods were dupli- 
cated on the soybean graphics, and the limits of the critical soybean periods 
were duplicated on the corn graphics in order to depict this overlap. 

Tables 3-1 through 3-4 were compiled to show the relative times of occurrence 
of peak greenness and to determine the stability of the reference growth 
stages to peak greenness in the years 1978 and 1979. On the average, corn 
peaked 29 days earlier than soybeans. In view of the time duration of the 
critical maximim greenness periods of 30 days for corn and 25 days for 
soybeans, there would be an average overlap of 12.5 days between the periods 
for the two crops if the midpoint of each period were placed at exactly peak 
greenness. There were other factors, however, related to separation (or 
discrete identification) of corn and soybeans that needed to be taken into 
account for establishing relative locations of labeling windows (heretofore 
known as critical periods); these factors are discussed In section 4.4. This 
relative closeness of peak greenness showed, nevertheless, that a common logic 
could be used for identifying corn and soybeans as summer crops In the U.S. 
Corn Beit states and Minnesota. Farther south (at least in Mississippi), the 
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TABLE 3-1.- RELATIONSHIPS OF CORN PEAK GREENNESS TO STANDARD GROWTH STAGES FOR YEAR 1978 
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TABLE 3-2.- RELATIONSHIPS OF CORN PEAK GREENNESS TO STANDARD GROWTH STAGES FOR YEAR 1979 
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TABLE 3-3.- RELATIONSHIPS OF SOYBEAN PEAK GREENNESS TO STANDARD GROWTH STAGES FOR YEAR 1978 


ORIGINAL PAGE IS 
OF POOR QUALITY 



3-11 



TABLE 3-4.- RELATIONSHIPS OF SOYBEAN PEAK GREENNESS TO STANDARD GROWTH STAGES FOR YEAR 1979 
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crop calendars Indicated a much wider separation In peak greenness, which 
would require execution of separate logic steps for the two crops where 
significant acreages are grown. 

3. 3. 4. 3 Relationships of Critical Periods to Standard Growth Stages 

Tables 3-1 and 3-2 showed that there was relative stability during the two 
crop years between the times of occurrence of corn tassel ing and denting as 
related to peak greenness. Fifty-percent tassel Ing occurred each year about 
5 days earlier than peak greenness, and denting occurred about 35 days later 
than peak greenness. Crop calendar data for the two crop years also showed 
that 50 percent tassel Ing occurred approximately 41 percent of the way through 
the growing season; that denting occurred about 66 percent of the way through; 
and that maturity occurred about 77 percent of the way through. Provided that 
the length of the growing season is known, additional knowledge of either the 
planting date or the peak greenness date (determined spectrally) could enable 
one to make good approximations of the time of occurrence of the growth 
stages. 

Tables 3-3 and 3-4 showed a relative existence of stability during the two 
crop years between the times of occurrence of soybean podding and soybean peak 
greenness. Fifty-percent podding occurred approximately 10 days before peak 
greenness. The relative time of occurrence of turning was not as stable. 

Some other relationships observed were that podding occurred about 50 percent 
of the way through the growing season and that maturity occurred about 
87 percent of the way through. Using these relationships, some valuable 
approximations regarding occurrence of soybean growth stages can be made, 
provided that the length of the growing season Is known and either the 
planting date or the peak greenness date Is known. 

Based on the known dates for peak greenness for corn and soybeans and the 
seasonal tassel Ing, podding, and maturity relationships described In the last 
two paragraphs, the midpoints and time spans of critical periods were 
calculated for those critical periods not having acquisitions near the period 
extremes (and allowing interpretative approximations). Calculations of the 
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postharvest period were made for all segments, based on maturity dates, and 
confirmed by Interpretation, where possible. Real dates determined either by 
ground observation (planting and harvest) or by nominal crop calendar are 
unbracketed In tables 3-1 through 3-4. Those dates determined by calculations 
are bracketed. Fifty-percent planted dates In the tables were based on 
special fields observations, where possible. 

3.4 LOGIC FOR DISCRETE IDENTIFICATION OF CORN AND SOYBEANS 

The ground-truth-labeled scatter plots were examined for all dates before, 
within, and after the periods Identified In section 3.3.3 for maximum green- 
ness of the two crops, to determine whether separation of corn and soybeans 
(when they occurred together) could be accomplished using the same acquisition 
used to Identify the two crops as summer crops. The separability of the crops 
In each acquisition examined Is recorded on the graphics In appendix C. The 
separability of other summer crops such as alfalfa, sunflowers, and cotton was 
also studied. Green number versus brightness scatter plots of the acquisition 
best illustrating separation are shown In appendix D for the various suimier 
crop combinations. The following paragraphs explain the processes necessary 
to achieve separation (Identification) of corn and soybeans In each case. 

3.4.1 CORN AND SOYBEANS IN THE U.S. CORN BELT STATES 

Acquisitions that fall in the maximum greenness period for Identifying 
soybeans as summer crops serve best for separating corn from soybeans. 

However, an acquisition in which corn has senesced sufficiently to be past the 
bright red stage Is best. If the peak greenness points of the corn and 
soybean crops are too close, this condition will not occur before soybeans 
begin to senesce, and overlap between the spectral distributions will then 
prohibit separation (fig. 3-1). If the correct acquisition Is available (and 
selected), the bright red colors of soybeans and the dark reds and browns of 
the corn will be easily distinguishable In the PFC product 1; and a line 
separating the two crop distributions (colors) can be placed In the scatter 
plot for the acquisition. 


Green number Green number 


ORIGINAL PAuc IS 
OF POOR QUALITY 



(a) Ideal separation acquisition. 



(b) Corn peak too near soybean peak 
(or soybeans not yet peaked). 


Figure 3-1.- Good separation versus overlapping distributions. 
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Figures 3-2 and 3-3 further Illustrate how greenness and brightness combine to 
give Ideal separation. Note that greenness Is the predominant separator early 
In the separation period (window) but that brightness may predominate late In 
the window. Dots plotted are typical pixels taken from segment 0943 for 1978. 

As stated in the Introduction, only pure dots could be labeled for training 
the classifier. While pure dots were defined as being those that were 
entirely within a field on all acquisitions used In labeling, some additional 
pixels In the best separation acquisitions looked like Impure pixels or were 
atypical of corn and soybeans. These pixels were always lying outside the 
main dot clusters In the plots. After examining several plots, It was 
concluded that a width of 13 brightness counts spanning a line lying along the 
green arm would include all of the main cluster pixels and would exclude less 
than 5 percent of the tctal corn and soybean dots that were otherwise pure. 

It was therefore decided to put limiters In the plots to eliminate these 
atypical pixels from consideration In labeling, In the hope that this would 
enhance classifier training. 

It was also noted that most of the corn dots that would be labeled "N" by the 
green number threshold of 13 were also below the limits of the main corn 
cluster and thit the* looked very much like trees in the PFC product 1; this 
confirmed that the iarlier decision to use 13 as the threshold between crop 
and noncrop wts appropriate. 

3.4.2 CORN, CORN AND HAY, OR CORN AND SUNFLOWERS 

These combinations of corn and other crops can be present in those acquisi- 
tions which best Identify corn as a summer crop. Winter or spring grains may 
also be grown in areas bordering the U.S. Corn Belt where these combinations 
occur; but, by using logic already described, these latter crops should sepa- 
rate readily as being nonsummer crops. In areas where soybeans are not grown 
or in areas where the time separation between peak greenness of corn Is too 
great to allow identification of corn (as a summer crop) using a green-up 
acquisition for soybeans, it is necessary to use the optimal green-up 
acquisition for corn (selected in accordance with section 3.3.1). Alfalfa 
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Brightness 


f 



Time, months 

(a) Greenness through time. 



F MAMJJASOND 
Time, months 


(b) Brightness through time. 

Figure 3-2.- Corn and soybean separation window duration. 
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Figure 3-3.- Corn and soybean separation window duration illustrated by trajectories through ti*r; 


hay, which Is the l-.v most likely to be confused as a summer crop, will be 
much greener ."nd *r y ter than corn. This will also be true of sunflowers. 

In these scatter p. . or In plots where there Is corn alone, the color of 
corn Is dark red, and a boundary can be placed above the corn In the plots, 
separating Its cluster from any confusion. 

3.4.3 CORN AND SPRING GRAINS 

Another possibility 1;: *nd north of the II, S. Corn Pelt states Is for corn and 
spring grains to be grown without soybeans being present. When this Is true, 
the optimal green-up acquisitions for Identifying corn may also contain spring 
small grains with green numbers higher than the corn threshold of 13; this 
will Identify the spring small grains as being a summer crop. However, It was 
shown using segment 0185 In Minnesota that In acquisitions late In the 
critical corn period the green numbers of spring small grains tend to drop 
below 13. The logic "Is the green nir.ber greater than or equal to 13 on all 
acquisitions?" will retain most corn as a summer crop but eliminate most of 
the spring small grains, provided that one early acquisition and one late 
acquisition are picked within the critical period. The corn can then be 
Isolated In the scatter plot of the late acquisition In the same way as in 
section 3.4.2. 

South o the U.S. Com Belt, two acquisitions should not be required because 
spring small grains In more southerly states are earlier maturing varieties. 

3.4.4 SOYBEANS OR SOYBEANS AND SPRING GRAINS 

Segments have not been found In the U.S. Corn Belt states where soybeans are 
present without corn, although Missouri has some segments where the percentage 
of corn Is low (e.g., seqment 0209 has about 12.6 percent corn). In the 
critical maximum green acquisitions for soybeans (as selected In sections 3.3.2 
and 3.4.1), the spring small grains will be below the summer crop green number 
threshold of 13. The soybean cluster will be easy to Isolate by viewing the 
bright red color of soybean dots In the PFC product 1 and placing a line In 
the summer crop scatter plot just below their lower limit. 
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3.4.5 CORN, SOYBEANS, SPRING GRAINS, AND SUNFLOWERS 

Landsat segments containing corn, soybeans, spring small grains, and 
sunflowers are not known to exist south of latitude 43° N. In the United 

States. One segment was found north of latitude 43° N. In Minnesota 

(segment 0185) that had sufficient acquisitions at the right times to 
illustrate the separation of corn and soybeans when this combination of crops 
is present. 

As shown in section 3.4.3, use of both an early acquisition and a late 
acquisition in the maximum greenness period for corn will separate corn from 
spring small grains. This is accomplished by applying the logic "is the green 
number greater than or equal to 13 on all acquisitions?" to the two acquisi- 
tions. The logic holds true for corn, but not for the soybeans because they 
will not be sufficiently emerged in these acquisitions. 

In order to retrieve the soybean dots after identifying corn, it is necessary 

to use an acquisition late in the soybean maximum greenness time period and 
generate another scatter plot. All dots not previously called corn are 
subjected to the logic "is the green number greater than or equal to 13?", and 
only the dots meeting this criterion are plotted on the scatter plot. 

While the corn is gone, sunflowers as well as soybeans are now identified as 
summer crops. Inspection of the plot shows, however, that the sunflowers 
occupy the riormal position of corn on the plot and are separated (in greenness 
and brightness) from the soybeans. The soybeans will be bright red in the PFC 
product 1, and the sunflowers will be dark red. Placement of a line in the 
scatter plot immediately below the last dot that is bright red in the scatter 
plot correctly separates 100 percent of the soybeans with no sunflower 
confusion remaining. 

Although requiring three acquisitions in the maximum greenness periods for 
corn and soybeans and application of summer crop logic twice, separation of 
corn and soybeans in this rather confusing mix of spring crops and sunflowers 
is possible. This is an example in which usage of temporal information on 
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both greenness and brightness is shown to be an effective method for labeling 
crops that peak in greenness at about the same time. 

3.4.6 SOYBEANS AND COTTON 

The only available blind sites growing soybeans but not corn were located in 
Mississippi. Application of the same logic, “Is the green number greater than 
or equal to 13?", to the maximum greenness acquisition for soybeans resulted 
in effective identification of soybeans as a summer crop In the Mississippi 
segments (numbers 0195 and 0200). However, cotton also was called a summer 
crop. Trial and error placement of limiters in the scatter plot showed that 
80 percent of soybeans would be called summer crops, committing only two 
cotton dots in segments 0195 and 0200 (1978), if the distance between the 
limiters was widened to 15 counts. Only 69 percent was called summer crops 
using a lesser width of 13 counts. The necessity for the increased width 
between the limiters could be due to brighter soil or to existence of a 
different sun angle. Soils south of the U.S. Corn Belt generally contain more 
sand (ref. 7). 

It was noted that soybeans in Mississippi were not as bright red in PFC 
product 1 as they were in the U.S. Corn Belt, and otner more northerly states. 
It was not possible place a line perpendicular to the green arm in these 
segments that would clearly separate soybeans as a separate distribution along 
the green arm. The amount of corn, however, was insignificant (<2 percent) 
and had senesced to the point that it was no longer identified as a summer 
crop using the acquisition suited to soybeans. Therefore, there was no 
confusion with corn. 

3.4.7 HAZE PRECAUTION 

Available hazy acquisitions in the development data set that fell in separa- 
tion windows for any opticn were examined to determine the effect of haze on 
the amount of brightness spread in the crop clusters. 

It was evident that in some separation acquisitions haze caused some crop dots 
to be outside the boundary of the second limiter because of Increased 
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brightness width of the crop clusters. The plot for acquisition 230, 
segment 0828, is an example of this. As a general rule, if more than 
10 percent of the pure dots for either corn or soybeans are lost to labeling 
in this manner, the acquisition is too hazy for use. 

3.5 USE OF LABELING WINDOWS FOR DESIGNATING CRITICAL ACQUISITIONS 

Section 3.3 showed that the green numbers could be used in the periods of 
preemergence, maximum greenness, and harvest and postharvest to identify corn 
and soybeans as summer crops. There is sufficient overlap in some areas 
between the periods for the two crops to allow common acquisitions to be used 
for this purpose. Additionally, section 3.4 showed that acquisitions in the 
maximum greenness period can be used to make discrete identification of cor' 
and soybeans (after they are known to be summer crops) using green number 
versus brightness scatter plots. The times of occurrence of the critical 
periods were shown in section 3.3.4 to be rather stable when related to crop 
growth stages reported in historical crop calendars. It was decided, because 
of this year-to-year stability, to establish labeling windows at the times of 
the critical periods. 


The establishment of labeling windows would provide a simplified terminology 
for designating the time periods in which acquisitions must occur. Based on 
the study of the dots, the following approximate window relationships should 
be established. (Note: C = corn; Y = soybeans.) 


Corn 

Soybeans 
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The alphanumeric codes appearing in this diagram represent the following 
windows and subwindows: 

Planting windows: Cl and Y1 

Green-up windows: C2 and Y2 (subwindows C2S1, C2S2, Y2S1, and Y2S2) 

Postharvest windows: C3 and Y3 

The diagram illustrates that the planting and postharvest windows for corn and 
soybeans almost coincide and that the green-up window for soybeans is later 
than the corresponding one for corn. These windows are based on the critical 
time periods plotted in the appendix C graphics. Appendix E defines each 
window in terms of a specific number of days. The amount of overlap between 
the C2 and Y2 windows will vary according to location. Two independent sets 
of windows were established in order to increase flexibility of their use, 
depending upon geographic location and the crop mix. Subwindows were 
established in order to facilitate designation of early or late acquisitions 
within windows. 

Figure 3-4 illustrates the location of these windows in relation to greenness 
through time profiles of typical crops in and near the U.S. Corn Belt states. 
Note that window Cl/Yl greenness threshold of 10 separates summer crops from 
pasture and winter wheat (or winter grains). Window C2/Y2 threshold of 13 
identifies corn and soybeans as summer crops. Window C3/Y3 is seen as not 
being very critical, but occasionally it is effective in separating out trees 
not called "N" by the other logic. Cotton, the only other summer crop 
addressed but not shown in the figur?, peaks in greenness a little earlier 
than soybeans and separates using window Y2 brightness as described in 
section 3.4.6. 

Based on the window definitions and their possible applications, as explained 
In sections 3.3 and 3.4, it was now possible to establish a series of labeling 
options. Table 3-5 lists the possible crop mixes of major crops observed in 
the development segments and establishes an appropriate labeling option for 
each mix. The option designates particular logic and combinations of windows 
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TABLE 3-5.- OPTION SELECTION TABLE 


Major crops 
present 

( >5% of pseudocounty) 

General area of the United States 

Between 1 at. 
37° N. and 
43° N. 

North of lat. 
43° N. 

South of lat. 
37° N. 

Corn and soybeans 

1 

1 


Corn, soybeans, and 

1 

1 


spring grains 




Corn 

2 

2 


Corn and hay 

2 



Corn and sunflowers 


2 


Corn and spring grains 

2S 

2S 

2 

Soybeans 

3 

3 

5 

Soybeans and 

3 

3 


spring grains 




Corn, soybeans, anc 

3S 

3S 


sunflowers 




Corn, soybeans, spring 


3$ 


grains, and sunflowers 




Soybeans and 


4S 


sunflowers 




Soybeans and cotton 



5 

Summer crop (general) 

0 

0 

0 

Winter grains 

Included in 

included in 

Included in 


all options 

all options 

all options 
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greenness curves. 


to be used in each case. Options that use subwindows are designated with an 
"S" and are called split options. Blank spaces In the table Indicate that 
either (1) the crop combination listed Is not known to be grown at the 
particular latitude or (2) no segments were available to confirm the 
effectiveness of the option. For example, an option for processing corn and 
sunflowers south of latitude 37° N. was not provided because this combination 
of crops was not observed In any developmental segments at these latitudes; 
and no segments were available for corn south of latitude 37° N., although It 
was probably processable using option 2. Further study would be required to 
establish option effectiveness for the blank spaces. 

The latitude dividing lines of 37° and 43° N. were established to roughly 
enclose the states comprising the U.S. Corn Belt. It was also noted that, 
partly because of soil factors and crop types and varieties being grown, crops 
north and south of this latitude envelope needed slightly different labeling 
rules (options) to separate them from corn or soybeans when compared to rules 
for the U.S. Corn Belt. 

Options are designed strictly for labeling corn and/or soybeans, with or 
without the presence of other crops. Although some options might prove useful 
as tools in labeling other crops, further testing Is advised prior to making 
such a usage (e.g., for sunflowers). 

Because of sampling techniques used in the United States, which place segments 
primarily on cropland (pseudocounty), a lower limit of 5 percent was used as a 
cutoff for categorizing crops In the major (or significant) category. In an 
area where sampling is not as precise or In areas outside the United States 
where statistics are not as accurate (historical for region, province, etc.), 
a lower limit of 2 percent Is recommended because large fluctuations may 
occur. Such fluctuations may cause acreages In Individual segments to be 
quite large. 

Table 3-6 specifies the number ov acquisitions to be used In each window when 
applying e specific labeling option. Footnotes to the table explain that 
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TABLE 3-6.- WINDOW ACQUISITION* REQUIREMENTS FOR EACH LABELING OPTION 



‘consecutive days trill not be counted In determining minimum numbers of acquisitions or for 
establishing aqy maximum limitations. 

b The overlap between windows C2 and Y2 mst be <15 days, and the gap between these windows xust be 
<5 days In order for option 1 or 1C to be a valid option. If gap >6 days, option IT may be used 
but 1C may not. 

C A minimum of one acquisition is required In window Cl or Y1 In all options. In areas where 
potential com and/or soybean fields green up prior to spring plowing, a minimum of two acqul* 
sltlons Is required In order to Increase the probability of detecting the planting signature. 

A minimum of two acquisitions Is required In all com and soybean states except Iowa, Illinois, 
and Nebraska. 

^Optional; may replace an acquisition In another window or subwindow under certain conditions 
(to be exercised by analyst only; see expanded option definitions In appendix F). 

e 0ne or two acquisitions are desirable but are not required. 
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consecutive-day acquisitions are not to be used when totaling the acquisition 
requirements. Footnotes also key In certain restrictions alreacly explained In 
an earlier section. This table capsullzes the Information In appendix F, 
which additionally provides alternative option selections (for some options). 

Options that use corn windows are for areas that grow corn but no soybeans, or 
areas where there Is low overlap between corn and soybeans. 

3.6 LABELING OPTIONS AND DEVELOPMENT SEGMENTS; ACCURACIES ATTAINED IN THE 
MANUAL MODE 

For each labeling option previously defined, logic diagrams were constructed 
(appendix G). Using the green numbers In the necessary acquisitions, each 
diagram shows the steps necessary to arrive at Identification of summer crops 
for the option. Note that a single diagram m^y apply to more than one option. 
Summer crop dots (Z) from these options are then Identified as either corn or 
soybeans by applying the appropriate scatter plot technique from appendix D. 

An option was selected for each development segment, and the summer crop logic 
was executed manually for each pure dot In the segment. Table 3-7 shows the 
results obtained. Separation accuracy of the dots was then evaluated using 
the acquisitions giving the best separation and placing the decision boundary 
as described in section 3.4. See table 3-8 for these results. 

Note that accuracies of both summer crop Identification and final corn and 
soybean separation for 1978 are less than for 1979. The reason for this Is 
unknown, but the procedure may not compensate for such factors as wide 
variations In moisture or episodic events. The accuracies attained outside 
the U.S. Corn Belt compared well with accuracies In the U.S. Corn Belt states 
(the Missouri segment Is In a marginal location). The least accurate result 
was obtained for the segment In the spring small grains area (73 percent), and 
the majority of the error was caused by low com accuracy. However, there 
were only seven pure corn and eight pure soybean dots to be labeled In this 
segment. 
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TABLE 3-7.- LABELING ACCURACY OF THOSE PURE DOTS CONSIDERED 



Spring small grains 
Minn. 

Pure corn 
Nebr. 

Pure soybeans 
Mi ss. 


U.S. Corn Belt 


Pure soybeans 


Seg. no. 
(year) 


0127 (1978) 
0843 (1978) 

0809 (1978) 

0141 (1978) 
0893 (1978) 

0205 (1978) 


Average 


0185 (1978) 


0222 (1978) 


0195 (1978) 
0200 (1978) 


Average 


0127 (1979) 
0209 (1979) 
0828 (1979) 
0893 (1979) 


Average 


0195 (1979) 


Z:Z, a 

% 
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TABLE 3-8.- SEPARATION OF CROPS CALLED SUMMER CROPS 


(a) 1978 separation of crops within limiters 


Area 

Seg. no. 

Corn, % 

Soybeans, % 

U.S. Corn Belt 

0127 

100 

95.5 


0843 

100 

100 



93.5 

87.5 


0141 

100 

100 


0893 

100 

96 


0205 

89 

60 

Average 

97 

89 

Spring small grains 

0185 

100 

87.8 

1978 Average 

97 

89 


(b) 1979 separation of crops within liniters 


Area 

Seg. no. 

Corn, % 

Soybeans, % 

U.S. Corn Belt 

m 

98.4 

100 



100 

100 


IS 

100 

100 


0893 

100 

100 

1979 Average 

100 

100 
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3.7 PREDICTION OF TIME OF OCCURRENCE OF LABELING WINDOWS; ACQUISITION 

SELECTION 

The goal of acquisition selection has always been to select those acquisitions 
which give the best spectral separation of the crop of Interest. For labeling 
purposes, this translates Into those acquisitions having the necessary 
characteristics for crop Identification. Features recognizable In the Imagery 
using Imagery Interpretation methods may be "gauged" using green numbers or 
other indices of vegetation In order to make this Identification using the key 
acquisitions. Some methods available for selecting these acquisitions, with 
pros and cons, are as follows. 

3.7.1 HISTORICAL CROP CALENDARS ADJUSTED BY INTERPRETATION 

In this method, growth stages are assumed to relate directly to labeling 
windows. Historical crop calendars, showing key growth stages on one or more 
years of data, are compared to the PFC Imagery products to determine whether 
the crop of Interest appears to be at the correct historical growth stage on 
the date of a given Landsat Image. The Interpreter uses his experience (using 
one or all available acquisitions) to Identify the crop of interest and to 
interpret the approximate growth stage of the crop of Interest in each 
acquisition. Selections of the acquisitions most nearly having the necessary 
crop characteristics are then made. 

This method, although giving the analyst the satisfaction of feeling that he 
is exerting direct Influence on the outcome of labeling and classification, i> 
very prone to error. Reference 8, which addresses some of the weaknesses of 
the procedure used during the Transition Year, attributes much of the labeling 
error to using this method of acquisition selection. The main weakness Is the 
average analyst's inability to recognize all of the growth stage variability 
in the crops. Also, inaccuracies In the historical crop calendars often 
confuse the analyst. 


3-31 


3.7.2 HISTORICAL AVERAGES OF GROWTH STAGE OCCURRENCE 


This method uses mc-an historical averages of growth stage occurrence as the 
basis for window establishment. A statistic for past years, such as 
50-percent planted. Is taken as the reference for locating a window. fhls 
statistic may be taken either from a crop calendar or from reliable 
statistical tables. This approach works well for years In which the crop of 
Interest Is planted In a normal time frame and no abnormal events occur 
affecting the growth cycle of the crop. 

A small test of the probable success of using this approach for 1980 was made. 
The average date of occurrence of midpoints of the two most critical crop 
windows, C1/C2 and Y1/Y2, wss calculated for segment 0893 corn and soybeans, 
using 1978 and 1979 crop calendars. The a erage dates obtained and compared 
to dates taken from 1980 calendars were as follows: 

Midpoint 1980 

Window average date calendar dat e Pi fference 


Cl 

135 

118 

-17 days 

C2 

196 

188 

-8 days 

Y1 

139 

136 

-3 days 

Y2 

22 1 

224 

+3 days 


In option 1 segments, such as segment 0893 where the windows for both crops 
are based on the latest crop planted (l.e., soybeans), the statistical average 
method might work well. However, In an option 2 segment, windows based on 
this example might have resulted In omission of a considerable amount of corn. 
If segment 0893 had a window Cl midpoint of 135 for 1980, this would mean the 
limits cf the window would range from day 115 to day 156. Tf some corn were 
planted as early as day 115 and the real window ranged from day 098 to 138, as 
Indicated by a midpoint of 118 on the 1980 crop calendar, an acquisition 
acquired between days 139 and 156 could result in much of the corn that had 
already emerged being thrown out as natural vegetation (using the established 
green number threshold of 10). The window for planting would, in fact, be 
overlapping onto emerged corn. 
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3.7.3 USE OF THERMAL AND PHOTOTHERMAL PHENOLOGICAL MODELS 


The Cross-Zuber thermal model for prrdlctlng corn growth stages (ref. 9) and 
the Majors-Johnson photothermal model for predicting soybean growth stages 
(ref. 10) were scheduled for evaluation as part of this developmental task for 
labeling corn and soybeans. If It was possible to use or adapt the growth 
stages from these models for selecting acquisitions, it was desired to do so. 

Certain nominal, year-specific crop calendar growth stages were shown in 
section 3.3.4 to be rather stable from year to year with respect to peak 
greenness. Windows calculated using peak greenness adequately encompassed the 
necessary acquisitions. Taking window midpoints based on these windows as the 
reference for determining the accuracy of stages determined using the 
Cross-Zuber and Majors-Johnson models will, in effect, give an indication of 
the stability of the model with respect to peak greenness (i.e., peak 
greenness becomes a common denominator for growth stage determination, window 
limits, and model evaluation). 

All of the developmental segments for which weather data were readily 
available were run using the two models. Tables 3-9 and 3-10 show the results 
of the Cross-Zuber model predictions and, also, the resulting bias when 
comparing the predictions to the stable crop-calendar-derived windows. The 
Cross-Zuber stages used as the basis for window midpoint prediction were: 

Window Cl = 50 percent planted (+1 day) 

Window C2 = 50 percent silked (comparable to 50 percent tassel ing) (+5 days) 
Window C3 = 50 percent mature (+45 days) 

The mean error in predicting the two most critical windows. Cl and C2, for 
1978 was 7 to 8 days and for 1979 it was 4 to 8 days. For window C3, the mean 
error for 1978 was 13 days and for 1979 it was 8 days. On the average, the 
model was early for 1978 and late for 1979. 

Tables 3-10 and 3-11 show the results of running all of the segments with the 
Majors-Johnson prediction model, for which weather data were on line. Also 
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TABLE 3-9.- COMPARISON OF CROSS-ZUBER MODELED GROWTH STAGES 
TO MANUALLY DETERMINED WINDOW MIDPOINTS FOR YEAR 1978 


State 

Seg. 

Window Cl 

Window C2 

Window C3 


no. 

Man. 

Mod. 

Bias 

Man. 

Mod. 

Bias 

Man. 


Bias 

Ind. 

0127 



FI 



-1C 

HI 

292 

-18 


0843 


140 

Sfl 

■ 

w 

-7 

m 

300 

-15 

in. 

0809 

135 

138 

+3 

205 

210 

+5 

309 

301 

-8 

Iowa 

0141 

156 

135 

wm 

214 

205 

-9 

307 

291 

-16 


0893 

138 

133 

H 

213 

202 

-11 

302 

288 

-14 

Mo. 

0205 

138 

139 

H 

217 

206 

-11 

302 

293 

-9 
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TABLE 3-10.- COMPARISON OF CROSS-ZUBER MODELED GROWTH STAGES 
TO MANUALLY DETERMINED WINDOW MIDPOINTS FOR YEAR 1979 


State 

Seg. 

Window Cl 

Window C2 

Window C3 


no. 

Man. 

Mod. 

Bias 

Man. 

Mod. 

Bias 

Man. 

Mod. 

Bias 

Ind. 

0127 

131 

138 

+8 

209 

214 

n 

300 

311 


Mo. 

0209 

135 

140 

+5 

218 

214 

H 

311 

305 

-6 

111. 

0828 

135 

138 

+3 

199 

212 


298 

305 

+7 

Iowa 

0893 

134 

143 

-1 

202 

218 


303 

320 

+17 




.75 


.5 

Mean: +7 

.25 




.0 


.8 

RMSE: 11 

.1 
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TABLE 3-11.- COMPARISON OF MAJORS-JOHNSON MODELED GROWTH STAGES 
TO MANUALLY DETERMINED WINDOW MIDPOINTS FOR YEAR 1978 


State 

Seg. 

Window Y1 

Window Y2 

Window Y3 


no. 

Man. 

Mod. 

Bias 

Man. 

Mod. 

Bi as 

Man. 

Mod. 

Bias 

Ind. 

m 

E9 

wm 

-1 

239 

231 

-8 

319 

293 

-26 


Hi 

Si 

KgV 

-1 

240 

233 

-7 

298 

294 

-25 

111. 

0809 

134 

147 

-13 

222 

232 


298 

294 

-4 

Iowa 

0141 

149 

eg 

-4 

234 

228 

-6 

316 

290 

-26 


0893 

149 

ESI 

-8 

228 

225 

-3 

316 

286 


Mo. 

0205 

140 

149 

+9 

230 

230 

0 

307 

293 

-14 



wmm 


ram 

.33 

Mean: -20.8 



BUflE 

Eli 


.55 

RMSE: 22.7 


TABLE 3-12.- COMPARISON OF MAJORS-JOHNSON MODELED GROWTH STAGES 
TO MANUALLY DETERMINED WINDOW MIDPOINTS FOR YEAR 1979 


State 

• 

Window Y1 

Window Y2 

Window Y3 


Man. 

Mod. 

Bias 

Man. 

Mod. 

Bias 

Man. 


Bias 

Ind. 

0127 

137 

150 

+13 

228 

234 

+6 

301 

299 

-1 

Mo. 

0209 

141 

151 

+10 

231 

234 

+3 

301 

300 

-1 

111. 

0828 

139 

148 

+9 

220 

233 

+13 

299 

297 

-2 

Iowa 

0893 

136 

151 

+15 

224 

236 



+12 

295 

304 

+9 



Mean: +11.8 

Mean: +8.5 

Mean: +1 

.25 



RMSE: 12.0 

RMSE: 9 

.5 

RMSE: 

4.7 
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shown is the resulting bias when comparing the results to the stable 
crop-calendar-derived windows. The Majors -Johnson model stages used as the 
basis for midpoint calculation were as follows: 

Window Yl = 50 percent planted 

Window Y2 * 50 percent pod beginning to fill (+12 days) 

Window Y3 = 50 percent mature (+38 days) 

The error in predicting the window Y1 and Y2 midpoints for 1978 was about 1 to 
2 days and for 1979 it was 8 to 12 days. For window Y3, the error was on the 
order of 1 and 21 days, respectively, for the 2 years. 

The mean error for windows 1 and 2, for both corn and soybeans, was not 
surprisingly high; but the root mean square error (RMSE) for six out of eight 
of the window determinations was above seven. This led to further examination 
of the error in the individual segments. Some bias in the individual segments 
was 2 weeks or more. Bias in three out of four of the window 3 determinations 
was also exceptionally high. 

The effect on labeling accuracy for two of the segments with the highest 
soybean window (Yl, Y2, and Y3) midpoint bias was determined. Segments 0809 
for 1978 and 0127 for 1979 were selected (both option 1). Comparison of the 
window limits to the data showed the following: 

Segment 0809 (1978): Window Yl limits of 122 through 163 (bias of 131 makes 

no selection of a window Yl acquisition and the segment 
remains unprocessable. Note: If an acquisition had 

been available on day 163, there would have been some 
omission because an omission of 12 percent corn and 
15 percent soybeans is indicated for day 164, which is 
1 day later than the cutoff. (See pages C-5 and C-ll.) 

Window Y2 limits of 231 through 255 (bias of 10) would 
cause the acquisition for day 244 to be selected for 
separating corn and soybeans. There is an additional 
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loss of corn of 16 percent with the threshold of 13, but 
no additional loss of soybeans when the acquisition is 
selected. Separation of the remaining corn is 
96.6 percent and of the remaining soybeans is 
94.1 percent. 

Window Y3 limits of 276 through 311 (bias of -4) would 
select no bad acquisitions. 

Segment 0127 (1979): Window Y1 limits of 125 through 165 (bias of 13) will 

cause this segment to become unprocessable, since two 
acquisitions In the window are not available. Day 138 
taken alone results in green number loss of about 
8 percent of the soybeans, but no corn is lost. 

Window Y2 limits of 233 through 257 (bias of 6) would 
result in selection of the day 247 acquisition. Green 
number loss of corn of 16 percent would result using 
this acquisition. Separation of the remaining crop dots 
was 100 percent. 

Window Y3 limits of 281 through 316 (bias of -1) would 
not result in the selection of any bad acquisitions. 

After evaluating these and other similar segment results, it was concluded 
that acquisitions selected using the models alone were unsatisfactory. 

When using the models, the amount of variation from segment to segment within 
a given crop year did not differ greatly from the variation observed when 
using the peak greenness (manual) method, as shown in table 3-13. 

The within-year variation of window midpoints in the U.S. Corn Belt states is 
sufficiently close that either a manual method, using maximum greenness and 
fractional season relationships, or the two phenological models should serve 
equally well if a seasonal correction could be applied to each window selected 
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TABLE 3-13.- AVERAGE DEVIATION FROM THE MEAN OF WINDOW MIDPOINTS 


Window 

1978 

1979 

Manual 

Model 

Manual Model 

Cl 

8 

3 

1 

2 

C2 

3 

2 

7 

2 

C3 

4 

4 

4 

5 

Y1 

6 

3 

2 

1 

Y2 

6 

2 

3 

1 

Y3 

7 

2 

3 

2 


by the models. The low variation shown in table 3-13 indicates that a single 
seasonal correction for all segments would be sufficient. 

3.7.4 USE OF PEAK GREENNESS TO CORRECT PHENOLOGICAL MODEL OUTPUT 

A method using peak greenness, while requiring no type of adjustment, would be 
limited to those segments having good acquisitions throughout the growing 
season. Numerous acquisitions are necessary to enable construction of 
accurate greenness curves for all of the selected dots used in the 
methodology. 

Segments can be labeled with fewer acquisitions than are needed to construct 
accurate greenness curves. In order to extend the accuracy of a peak 
greenness method of window determination to these additional segments, it was 
decided to see whether a seasonal correction applied to the output of the two 
phenolog’cal models could do this. This correction would be based on window 
midpoint dates of segments having acquisition histories sufficient to 
determine and apply peak greenness. As a test of this method, the sign on 
each mean bias determined In tables 3-9 through 3-12 was changed, and the 
quantity obtained was then applied as a correction to each model -determined 
midpoint in the appropriate table. (Remember that section 3.3.4 showed a 
direct relationship between crop calendar stages and peak greenness.) New 
window limits for each segment were then calculated, and the charts in 
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appendix C were checked to determine whether any undesirable acquisitions 
would be selected by using these new windows. No bad acquisitions resulted. 

Appendix H documents the details of a technique for determining windows using 
peak greenness and for applying the results of this technique In order to 
correct phenological models. 

3.7.5 CURRENT-YEAR CROP CALENDAR ADJUSTED BY INTERPRETATION 

This method Is similar to the method using historical crop calendars, and it 
has the same weaknesses when applied on a segment-by-segment basis. It has an 
additional weakness in that considerable resources are required to build 
current-year crop calendars in order to process the Landsat data for that 
year; this means that no segments can be processed until after these calendars 
have been constructed. 
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4. PROCEDURE ASSEMBLY 


4.1 MAN/MACHINE RELATIONSHIPS 

The procedure was supposed to be assembled in such a way that the objective of 
labeling a few pure pixels with high accuracy :ould be accomplished with the 
minimum number of manual operations and with the maximum utilization of the 
computer being made. Also, any analyst decisions to be made would be easy 
ones having objective decision rules, and they would be computer supported 
whenever possible. 

Assembly of the major procedural functions in modular form would make it 
easier to insert major Improvements at a later date; for example, replacement 
of a manual step with a fully computerized one. 

In order to make use of available computer programs, to establish the desired 
man and machine relationships, and to combine groups of manual and machine 
functions into modular form, figure 4-1 was constructed. Figure 4-1 depicts 
two machine-processing sessions that produce reports to be utilized by the 
analyst during two subsequent imagery contact sessions. These four sessions 
result in data and training labels, which are used in the final classification 
session. 

4.2 PROCESSING FUNCTIONAL FLOW 

The following subsections explain the purposes of each processing session, the 
rationale for use and integration of existing methodology, and the factors 
affecting the decision to create new methodology or software. 

4.2.1 SESSION 1: INITIAL ACQUISITION SELECTION 

This session, a computer processing step, uses the Cross-Zuber and Majors- 
Johnson model growth stage predictions as the basis for determining the 
optimal time periods for use in the pure-pixel isolation and profile parameter 
creation software and for establishing labeling windows. Model growth stages 
are used directly for determining the pure-pixel isolation period and the 
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Figure 4-1.- Functional flow chart. 
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profile parameter periods; a correction (appendix H) Is applied to the modeled 
growth stage In order to determine limits for the labeling windows. 

The program for running the phenol oglcal models, available on Laboratory for 
Applications of Remote Sensing System (LARSYS), was utilized by personnel of 
the Design and Integration Section to create a data base of growth stages 
based on data from the primary weather reporting stations. The data base Is 
now available for use on JSC computers and Includes the States of Iowa, 
Illinois, and Indiana. The software for running session 1, AIREPORT, accesses 
this data base and creates two reports: the Crop Report and the Corn and 

Soybean Window Selection Report. Examples of these products are contained In 
appendix I. 

The Crop Report shows the growth stages predicted by the two models. For the 
Corn and Soybean Window Selection Report, the corrections shown In table 4-1 
are applied to the appropriate model growth stages to arrive at procedure 
labeling window midpoints for 1978, 1979, and 1980: 

TABLE 4-1.- CORRECTIONS TO BE APPLIED TO MODEL GROWTH STAGES 



Corn 

Soybeans 

Year 


50% 

silked 

50% 

mature 

Bg 

50% 

podding 

50% 

mature 



+7 days 

+13 days 

-2 days 

mm 



-4 days 

-8 days 

-8 days 

-11 days 




-16 days 

-4 days 

+6 days 

-8 days 




a The corrections for 1980 growth stages are not based on the 
technique described in appendix H; but current-year crop 
calendars were used and adjusted by Imagery Interpretation 
(by the procedure developer), because of system problems 
In generating green numbers. 

Based on the corrected windows, the amounts of overlap and gap between windows 
are applied In order to determine which corn or soybean options can be 
applied; and based on the latitude of a segment, the report gives a complete 
listing of the options executable for a segment. 
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4.2.2 SESSION 2: IMAGE SCREENING, PROCESSABILITY DETERMINATION, AND 

OPTION SELECTION 

Error messages In the Corn and Soybean Window Selection Report will Indicate 
when segments are nonprocessable based on location, window overlap, or 
unscreened acquisition availability. If the segment Is still processable 
after session 1, the analyst proceeds with session 2. 

Images are screened In this session for usability. Based on the availability 
of usable acquisitions to fulfill requirements of later procedural sessions, a 
labeling option Is selected and certain key acquisitions are designated for 
use In running these later sessions. 

Based on the good labeling results obtained using the experimental logic 
(table 3-7), It was decided to maintain the 5-percent criterion for clouds, 
cloud shadows, or bad data. The analyst screens each of the acquisitions that 
are flagged by program AIREPORT (for potential use) to omit any acquisitions 
not meeting these data quality standards and to omit any acquisitions having 
misregistration of greater than two pixels. The acquisitions remaining after 
screening are Inventoried to determine whether the following requirements are 
met: 

• A minimum of two and a maximum of five acquisitions for use in program 
AUTOFLD, which selects pu*e field pixels. 

• A minimum of one preemerged and three emerged acquisitions, or a maximum of 
two preemerged and six emerged acquisitions, for use in programs AUTOCLS 
COEFGB and AUTOCLS COEFGBTR (session 4). These programs create a file of 
profile parameters that Is used In cropland and noncropland determination 
and In classification. 

• Acquisitions to meet the minimum labeling requirements determined using 
table 3-6. 

If the above acquisitions are available, the segment Is processable. A 
maximum of 12 acquisitions may be used because of spectral aid program 
limitations. Care must be taken to assure that there are no more than two 


4-4 


preemerged and six emerged acquisitions for use In AUTOCLS COEFGBTR. Because 
this program will automatically eliminate excessive acquisitions, key 
acquisitions used In labeling may be eliminated. It Is essential that as many 
labeling acquisitions as possible also be used in program AUTOCLS. because the 
labels are used to label those profile parameters obtained from AUTOCLS 
COEFGBTR for use In classification. Rules were devised whereby excess 
acquisitions that met the Image quality standards would be eliminated by the 
analyst rather than by the computer. 

In order to select a labeling option, It Is necessary to know the geographic 
latitude- of the segment and the major crops being grown In the area (see 
section 3.5). In the United States where historical data bases are readily 
available at the county level, It Is recommended that the data bases be used 
to determine crop percentages. When processing data on other countries, It 
may be necessary to use Information that Is In fragmentary and Incomplete 
form. In each case, the source and quality of the Information will need to be 
evaluated prior to Its use. The purpose In using this Information Is, of 
course, to enable the analyst to select the option that will best Identify 
corn and soybeans as summer crops and eliminate (label N) any confusion 
crops. High concentrations of crops that are not separable can cause highly 
erratic labeling and classification results. The best Information available 
should therefore be used in order to make the option selection. 

When a suitable cloud-screening algorithm can be designed and automated, and 
if temporal registration of acquisitions can be improved, much of this manual 
session can then be combined with session 1, thus Increasing efficiency. 

4.2.3 SESSION 3: SUMMER CROP AND OTHER LABEL CREATION; ANALYST 

AID PRODUCTION 

Acquisitions selected in session 2 and information on the worksheet from that 
session are used to create the necessary files for use In the session 3 
computer programs. 
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Program AUTOFLD (FLDPROC) Is run to Isolate pure pixels and create a file of 
dots (PURPIX file) to be labeled; program AUTOCLS C0EF6B creates greenness 
profiles through time and determines the a, a, and t Q parameters In the 
following equation and creates a file: 

G(t) • S 0 (t/t 0 )*e , ( t2 ' t o ) 

where 

G 0 = the soil line greenness 

t Q = the cdlculated spectral emergence date 

a, a = constants 

The FLDPROC program uses the edge gradient enhancement technique described In 
appendix J. This program and the AUTOCLS COEFGB program were adapted for 
multitemporal use In corn and soybean areas under the technical direction of 
Dr. G. Badhwar of NASA/JSC. Appendix M shows a flow chart of tnese programs. 

The remaining programs for session 3 process the pure pixels to determine 
which are summer crops, and then the summer crops labels are printed out on a 
spectral scatter plot and on green number versus brightness listings, as 
follows: 

1. Using the acquisitions indicated on the session 2 worksheet, the appropri- 
ate logic for the selected summer crop labeling option (appendix G) is 
executed in the CROP program. A file (LISTCROP) is created containing 
crop definitions (Z or N). 

2. The file of parameters from COEFGB Is then used in a program called PEAKGB 
to determine crop/noncrop by a profile parameter method (see appendix K) 
that uses a length of season parameter (o). Only the pixels in the PURPIX 
file that have survived the CROP program as summer crops are passed 
through this screen. The PEAKGB program is run in series to the CROP 
program as a greenness curve screen that does not depend on soil lines (as 
d r os the CROP logic with green numbers). Usually, CROP Is slightly more 
accurate; but in those cases where soil lines are inaccurate, the PEAKGB 
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program Is superior for detecting trees. A file (PURECROP) of summer crop 
(Z) and nonsummer crop (N) dots Is created for later editing (at the 
beginning of session 5). 

3. The summer crop labels from step 2 are also used as Input for the 

following programs that create spectral aids and listings. The labeling 
aids are output for use by the analyst In session 4. 

e SRCSTBLA - produces table A (example In appendix L, page L-2) 

e SRCSTBLB - produces table B (example In appendix L, page L-3) 

e SRCSTBLC - produces table C (example In appendix L, page L-4) 

e SRCSSPA - educes an unlabeled scatter plot of all 418 grid dots used 

In labeling dot selection or a scatter plot with only the pure summer 
crop dots (examples in appendix L, pages 1-5 and L-6) 

e SRC3TPA - creates green number through time and brightness through time 
trajectory plots (example In appendix L, page L-7) 

All of the programs to be executed In this session, except AUTOCLS COEFGB, 
were designed specifically to support this labeling procedure and the profile 
parameter classification system. The listings and plots, adapted from 
Procedure P1A, enable the analyst to quickly find a scatter plot dot In the 
Imagery, or vice versa, and to go from one listing to another for making 
various comparisons. Instructions for executing these programs either 
Individually or as combined runs are contained in reference 11. Shortcuts for 
the creation of the necessary files will be readily apparent to the computer- 
oriented analyst after becoming familiar with the profile parameter system. 

4.2.4 SESSION 4: CORN AND SOYBEAN LABEL ASSIGNMENT ■, APPROVAL CF N LABELS 

Prior to proceeding with final labeling. It is necessary at this point to 
check the accuracy of the soil lines because the green numbers based on them 
affect accuracy of the automated summer crop logic used In the programs for 
session 3. The same checks as those used In the development research (see 
section 2.4) were Inserted Into the procedure at this point, with Instructions 
for making soil line corrections. If necessary. If a reliable, error-free 
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I 

method for generating soil lines can be devised that will work In segments i 

having very small amounts of bare soil, this step may then be eliminated. 

The final (discrete) Identification of corn and soybean crops Is accomplished 
using the scatter plot of summer crop dots. Limiters are placed In this plot 
to identify the zones containing the crop dot clusters, and a decision 
boundary is placed between the corn and soybean cluster distributions. Since 
the labeled scatter plots may not always have sufficient uots within both the 

j 

corn and soybean distributions (e.g., option 2 or 5) to clearly define the 
slope of the green arm In the C2 or Y2 window acquisitions, rules for use of 
the unlabeled plot were devised to assist in establishing the first limiter 
(the second limiter is automatic upon establishment of the first one). The 
steps described in sections 3.4.1 through 3.4.6 were used for final placement 
of the line distinguishing the* corn and soybean distributions; a few simple 
rules were found to apply to multiple combinations of crops. 

Rules were devised for editing (checking purity) dots of all categories, C, Y, 
or N, prior to accepting them. In order to ensure that the proper spectral 
training would be provided for the classifier, C or Y dots needed to be on the 
same field in the window 1 and window 2 acquisitions. Since the value of the 
window 3 acquisition is found in the occasional labeling of an N dot, N dots 
had to be registered to the same field in windows 1, 2, and 3 in order to be 
accepted. 

After purity is checked, the analyst checks wnich crop distribution the dot is 
in to make final categorization of the dot. Between 15 and 20 dots from 
categories C, Y, and N are required for classifier training. 

A final action to be performed in this session is to select a pure blob of at 
least 10 ; e contiguous coi r. pixels to act as a starter for program COEFGBTR. 

A ulob should be picked and coordinates of Its vertices should be recorded in 
line-pixel format. 
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Various methods have been tried In order to automate the discrete Identifica- 
tion of corn and soybeans. However, the manual method of separating 
distributions described here has been shown to give the best results so far 
when using profile parameter classification. One means of further automation 
would be through the use of a cathode-ray tube display of the separation 
scattei plot. If dots In the displayed plot were colored the same color as In 
the PFC product 1, then cursorlng of a decision boundary could be made to be 
very fast and accurate. Software could then sort (and record) the dots 
according to crop category depending upon placement of the cursored line, and 
the resulting labels could be transferred to the classifier automatically. 

This method would give less accurate classification results, unless pure dots 
were always Isolated and plotted. (At present, purity Is manually confirmed.) 

4.2.5 SESSION 5: CLASSIFICATION (TWO OR THREE CATEGORIES) 

The labels accepted during session 4 are final edited by the execution of the 
PPIXEDIT program. The analyst enters the dot number and crop symbol for each 
of the pure labeis. The program transforms all dot synods into the required 
numeric labels for program PPIXGT. This results In a new file PURECROP. This 
file ma*y be printed If evaluation of the labels is desired at a later date. 

Before classification can be run, profile parameters for classification must 
be run using program AUTOCLS COEFGBTR. The flow chart in appendix M shows the 
flow of this program. The blob of starter corn pixels and other control cards 
are Input from file USERGBTR. Because of the length of time necessary to run 
this program, it Is desirable to begin running It as early as possible; and it 
can be run as soon as the procedure has identified the corn distribution In 
session 4 (In order to obtain a valid corn blob). 

Classification (program PARCLS) may be run after the necessary files arc 
available. The labels, profile parameters, some control cards, and symbol 
definitions are input. See appendix M, page M-6, for an illustrative flow 
chart. Details of the control cards are contained In reference 11. Either 
the three categories C, v , and N or a combination of any two of the categories 
may be classified, depending upon the availability of category dots. If a 
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category is determined to be insignificant (<7 dots), it may be disregarded. 

If there are between 7 and 14 of either C or Y dots, this category may be 
combined with category N. If both categories C end Y have between 7 and 14 
dots, the segment is nonprocessable. 

A report of classification and a classification map are printed as final 
products. The percentage of a crop in the scene may be determined by dividing 
the number of pixels classified in each category by 22,932. 
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5. TEST RESULTS 


5.1 GENERAL 

The procedure was tested in the semi automated mode to verify that satisfactory 
results could be obtained using an independent set of segments for the years 
1978 and 1979 In the U.S. Corn Belt states. For this test, the segments in 
table 2-2 were used. Personnel used for this test had the following 
qualifications: 

1. One Landsat analyst with no previous contact during procedure development 
for labeling 

2. One programmer used in procedure programming for computer processing 

3. The author (a Landsat analyst) for quality assurance 

The average time required to process a segment, executing sessions 1 through 
5, was 85 minutes analyst time, 70 minutes central processing unit (CPU) time, 
70 minutes operator time, and 22 minutes quality assurance l ~«e. Set 
figure 5-1. 

5.2 TEST RESULTS 

Appendix N contains an overall evaluation and individual evaluations for the 
seven segments. The aspects evaluated were classification accuracy, dot 
labeling accuracy (with general characterization of errors), and purity of 
labels sent to the classifier. The highlights of the results were as follows: 

1. Proportion estimation accuracy (area of seven segments) 

• Corn: 99.38 percent of ground truth (-0.62 percent relative error) 

• Soybeans: 95.57 percent of ground truth (-4.43 percent relative error) 

§ Other: 102.7 percent of ground truth (+2.73 percent relative error) 


/ 
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Figure 5-1.- Segment processing time 






















































2. Labeling accuracy (all dots In seven segments) 

• Corn called corn: 86 percent 

§ Soybeans called soybeans: 90 percent 

• Other called other: 98 percent 

3. Purity of labels after analyst acceptance: 95.8 percent 

Classification yielded relatively clean classification maps with an average 
percentage (of the area inventoried) of 68.45 percent of the pixels being 
correctly classified. 

Although a detailed characterization of classification errors was not made, 
the segments with higher errors seemed to show some correlation with labeling 
error causes. The majority of the labeling error was procedural ly inherent 
(i.e., error from factors that the labeling procedure, in the interest of 
simplicity, is not designed to handle). These errors appeared to be 
associated with moisture fluctuations, as evidenced by early and late fields 
or by episodic events such as hailstorms. 

5.3 LABELING ERROR CHARACTERIZATION 

Errors made in labeling the 152 corn (C), 143 soybean (Y), and 125 nonsummer 
crop (N) dots were as follows: 

1. C called Y = 10 dots 

5 - wrong distribution; possible late planting 

2 - incorrect decision boundary placement 

1 - clerical error 

2 - questionable ground truth 
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2. C called N = 11 dots 


2 - planting signature not detected 

5 - low green number (<13) in separation acquisition; probable early 
pi anting 

2 - sigma error 

1 - clerical error 

1 - questionable ground truth 

3. Y called C = 13 dots 

8 - wrong distribution; hail damage, moisture variations; possible poor 
stands 

2 - impure pixels 

3 - questionable ground truth 

4. Y called N = 2 dots 

1 - planting signature not detected 
1 - impure pixel 

5. N called Y = 2 dots 

1 - volunteer growth 
1 - mixed cropping 
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6. EXTENDABI L I TY 


Aside from such factors as soil differences, drought, pests, and plant disease 
effects that can affect light reflectance from crops, there appear to be some 
factors that could be considered common to all countries when adapting this 
procedure (or another) for Identifying crops and determining crop proportions 
using Landsat Imagery. Generally, these factors seem to be as follows: 

1. The crop has unique colors or geometry In the imagery products that permit 
visual (manual) identification of the crop, either In a data acquisition 
for a single date or In data acquisitions from multiple dates. All 
variations or varieties of the crop must be visually Identifiable. 

2. In order to develop automatic methods, some means (e.g., Kauth greenness, 
brightness) must be available for gauging the intensity of the identifying 
characteristic(s) in a manner that will identify the crop but exclude 
other crops or noncropland. It appears that automation accuracy can only 
be as good as manually executable accuracy (step 1), excluding clerical 
error. 

3. The time and duration of occurrence of the unique Identifying characteris- 
tics are relatable to the plants' growth cycle. 

4. A method must be available for accurately predicting the time of 
occurrence of the characteristics during each new crop year. 

5. Sufficient data acquisitions (at the right times) must be acquired in 
order to execute the desired logic for each segment. 

Because of the difficulty In identifying all variations of a crop, ground 
truth should be studied carefully when initially adapting the procedure and 
establishing decision rules. It was noted, for example, when developing the 
procedure, that soybeans in Mississippi did not look like soybeans in Iowa. 
Because of the greenness-brightness differences, part of the Iowa logic was 
executable In Mississippi, but not all of It. Since the differences were 
known, however, the procedure logic was adaptable. 


The relative times of occurrence of peak greenness were used In determining i 

the logic for identifying corn and soybeans In this procedure. If peak green- 
ness of a confusion crop occurs at a time very near that of the target crop, 
then greenness cannot be used to discretely identify the target crop unless ’ 

there are other unique features (e.g., brightness) that occur consistently and 
with sufficient, quantifiable Intensity. The two transformations of greenness 

w 

and brightness do have considerable flexibility of application in the multi- 
temporal sense. Adaptability of this procedure will depend upon relative 
times of greenness and brightness of various crops. 

Prediction of the time of occurrence of peak greenness (or other characteris- 
tics) using phenological models as ir. this procedure will depend largely upon 
availability of an appropriate means of starting the model at the right time, 
the amount of variability in planting date, or the variations of t.ie target 
crop. 


6-2 


7. CONCLUSIONS 


f 

1 
i 

« 
i 

The labeling accuracy of summer crops and the separation of corn and soybeans 
as Individual crops, when applying the logic to the research segments 
(section 3.6), Indicated a potential for procedural labeling accuracy slightly 
less than 90 percent for corn and soybeans for 1978, but slightly over 
90 percent for 1979. This, plus a potential of greater than 90 percent 
accuracy for nonsummer crops. Indicated that an assembled semlautomated 
procedure could provide the required training for a classifier. 

Testing of the assembled semi automated labeling procedure on seven Independent 
segments for 1978 and 1979 verified that sufficiently high classification 
accuracy (on the average) could be obtained using the labels. The average 
classification bias for the total area of the seven segments was -0.62 percent 
relative for corn and -4.43 percent relative for soybeans; Inputting labels 
with a labeling accuracy of 86 percent for corn, 90 percent for soybeans, and 
93 percent for other. 

The major part of both labeling and classification error seemed to be caused 
by precipitation variations and eplsodal events. Labeling error attributable 
to the semi automated procedure was approximately 6 percent, based on the test, 
and 4.3 percent was due to these two factors. 

The phenological growth stage models that were evaluated, with a spectrally 
determined seasonal correction, worked well at determining the correct 
labeling windows. The semi automated labeling logic, when applied to the 
window acquisitions, was then effective in arriving at correct labels. 

The requirement for high labeling accuracy slightly reduces the number of 
segments that can be classified because high labeling accuracy requires higher 
precision in acquisition selection. 

The more complexity there Is In numbers of crops grown, the greater the number 
of acquisitions and the more complex the green number logic required (see 
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section 4.4). Greater precision In acquisition selection (window width) Is 
also required. 

Provided that the following conditions can be met, the principles of this 
labeling procedure (and use of the profile parameter classification system) 
should be extendable to other countries: 

1. Stable time periods (windows) based on growth stages can be established 
for application of a semi automated green number logic to discretely 
Identify corn and soybeans. 

2. Additional logic can be devised for separating out any other summer crops 
(or confusion) and Identifying them as other. 

3. A phenological model is available for predicting the times of occurrence 
of the essential growth stages In new crop years. The model may or may 
not be seasonally adjusted, depending on the nature of the model. 
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APPENDIX A 


UNLABELED SPECTRAL AIDS 
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APPENDIX A 

UNLABELED SPECTRAL AIDS 

This appendix contains examples of unlabeled spectral aids used In procedure 
development. The following examples are given for segment 0843, located In 
Henry County, Indiana. 

1. First page of 209-dot green number listings for acquisition dates selected 

2. Unlabeled green number versus brightness scatter plot 

3. Example pages of greenness through tW and brightness through time 
spectral plots 
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LABELED SPECTRAL AIDS 

This appendix contains examples of labeled spectral aids used in procedure 
development. The following examples are given for segment 0843, located in 
Henry County, Indiana. 

1. First page of 209-dot green number listings for acquisition dates selected 

2. Scatter plots of planting (window Cl and Yl) acquisitions (dates 78151 and 
78160) with corn (C) and soybean (Y) dots plotted 

3. Scatter plots of window Cl and Yl acquisitions (dates 78151 and 78160) 
with other categories of dots plotted 

4. Scatter plot of green-up and separation (window Y2) acquisition 
(date 78232) with C and Y category dots plotted 

5. Scatter plot of window Y2 acquisition (date 78232) with other categories 
of dots plotted 

6. Scatter plots of postharvest (window C3 and Y3) acquisitions (dates 78305 
and 78313) with C and Y category dots plotted 

7. Scatter plots of window C3 and Y3 acquisitions (dates 78305 and 78313) 
with other categories of dots plotted 

The procedure processes the acquisitions in designated windows using a multi- 
temporal logic that separates the corn and soybeans from other categories of 
crops and noncropland. 
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APPENDIX C 


CRITICAL TIME PERIOD (WINDOW) CHARTS 

This appendix cental ns charts depleting the .'ata relevant to determining 
critical time periods (windows) for labeling corn and soybeans In the United 
States. The charts have he following captions: 

Figure C-l - Critical time period (window) chart for corn for year 1978 

Figure C-2 - Critical time period (window) chart for corn for year 1979 

Figure C-3 - Critical time period (window) chart for soybeans for year 1978 

Figure C-4 - Critical time period (window) chart for soybeans for year 1979 

Data recorded on each chart are as follows: 

1. For each segment, the segment number and location (county, state). 

2. The growth stages nearest the time of occurrence of the critical periods 
denned earlier. Dates of these stages are to be taken as points of 
reference for plotting the segment acquisitions on a relative scale. The 
"reference" stages are from the best available CRD or state crop calendar 
for the applicable year. Reference stages for the periods are: 

Crop Period Reference stage(s) 


Corn 

Preemergence 

50% planted 

Corn 

Maximum greenness 

50% tasseled, 50% dented 

Corn 

Harvest 

50% mature 

Soybeans 

Preemergence 

50% planted 

Soybeans 

Maximum greenness 

50% podded, 50% turned 

Soybeans 

Harvest 

50% mature 


3. Each acquisition date for each segment plotted in relation to the 
reference growth stages. 

4. The date of occurrence of growth stages recorded by the ASCS enumerator 
during the ground data surveys. Median dates were calculated whenever 
possible, but planting data were not always available. 
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5. Comments addressing the adequacy of acquisitions at early and late 
extremes of the critical periods. 

6. Approximations of the limits of the critical periods. These limits are 
empirically derived. The final limits recorded are based on calculations 
described in appendix H ( which are based on the relationships shown In 
tables 3-1 through 3-4. 

The following critical period symbology Is used in these figures: 

1. Dotted line brackets represent the corn window limits. 

2. Solid line brackets represent the soybean window limits. 

3. Dashed lines represent the corn/soybean window limits coinciding on same 
date. 
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Figure C - 2 .- Critical time period (window) cnart for corn for year 1979- 












ORIGINAL FAGS B 
OF POOR QUALITY 



Figure C-2.- Concluded 



ORIGINAL PAC r - K‘ 

» poor Quality 



C-9 


Figure C-3.- Critical time period (window) chart for soybeans for year 1978. 
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Figure C-3.- Continued 
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Figure C-3.- Concluded 
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Figure C-4.- Critical time period (window) chart for soybeans for year 1979. 
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Figure C-4.- Concluded 
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ILLUSTRATIVE CORN AND SOYBEAN IDENTIFICATION 
SCATTER PLOTS 

This appendix contains example scatter plots to Illustrate separation of corn 
and soybean distributions, when various combinations of summer crops are being 
grown. The crop key for crops grown in each segment Is recorded on each 
individual scatter plot. The following is a list, by figure nunber, of the 
various crop combinations shown: 

Figure D-l - Corn and soybeans In the U.S. Corn Belt states 
Figure D-2 - Corn or corn and hay 

Figure D-3 - Corn and spring grains or corn and sunflowers 
Figure D-4 - Soybeans or soybeans and spring grains 
Figure D-5 - Soybeans and sunflowers 
Figure D-6 - Soybeans and cotton 

In addition, an example of the effects of excessive haze on crop separation is 
shown in figure D-7. 
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Fiaure D-l.- Example scatter plot illustrating the separation of corn and soybean distributions 
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Figure D-3.- Example scatter plot Illustrating the separation of 
corn and spring grains or com and sunflowers (option 2S). 
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Figure D-5.- Example scatter plot ll^'itrating the separation 
of soybeans and sunflowers \ ^tion 3S or 4S). 




Figure D-6.- Example scatter plot illustrating the separation 
of soybeans and cotton (option 5). 
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CORN AND SOYBEAN LABELING WINDOWS 


Window Yl: 

Open: Soybeans 50 percent planted minus 25 days (Opening date may not be 

earlier than day 110.) 

Close: Soybeans 50 percent planted plus 16 days 

Window Y2: 

Open: Soybeans 50 percent podding minus 1 day 

Close: Soybeans 50 percent podding plus 23 days 


Window Y3: 

Open: Soybeans 50 percent mature plus 20 days 

Close: Soybeans 50 percent mature plus 55 days 

Window Cl: 

Open: Corn 50 percent planted minus 20 days (Opening date may not be 

earlier than day 110.) 

Close: Corn 50 percent planted plus 21 days 

Window C2: 

Open: Corn 50 percent tassel ing minus 10 days 

Close: Corn 50 percent tassel ing plus 20 days 

Window C3: 

Open: Corn 50 percent mature plus 30 days 

Close: Corn 50 percent mature plus 60 days 
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Subwindow C2S1: 


Open: Corn 50 percent tassel Ing minus 10 days 

Close: Corn 50 percent tassel Ing minus 1 day 

Subwindow C2S2: 

Open: Corn 50 percent tassel Ing plus 11 days 

Close: Corn 50 percent tassel Ing plus 20 days 

Subwindow Y2S1: 

Open: Soybeans 50 percent podding minus l day 

Close: Soybeans 50 percent podding plus 8 days 

Subwindow Y2S2: 

Open: Soybeans 50 percent podding plus 14 days 

Close: Soybeans 50 percent podding plus 23 days 
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APPENDIX F 


OPTION 0: 

Acquisitions 

Restriction: 

OPTION 1: 

Acquisitions 

Restrictions 


An acquisition 
consecutive day 


EXPANDED OPTION DEFINITIONS INCLUDING 
RESTRICTIONS AND ALTERNATIVES 

A. One window Y1 acquisition 1 in Iowa, Illinois, or 
Nebraska; at least two window Y1 acquisitions In other 
states. 

B. One window Y2 acquisition. 

C. One or two window Y3 acquisitions desirable but not 
requi red. 

Two-category labeling only: Z and N. Label 15 to 20 dots 

in each category. 


A. At least one window Y1 acquisition in Iowa, Illinois, or 
Nebraska; at least two window Y1 acquisitions in other 
states. 

B. One window Y2 acquisition. 

C. One or two window Y3 acquisitions desirable but not 
requi red. 

A. Overlap of windows Y2 and C2 must be <15 days. 

(Alternate option: zero.) 

B. Gap between windows Y2 and C2 must be <5 days. 

(Alternate option: If gap is 5 to 14 days, use 1Y.) 

C. Number of labelable corn dots must be >15. 

(Alternate option: 1Y.) 


fill be counted as nonconsecuti ve-day coverage only. Only one 
may be counted in these minimum date determinations. 
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OPTION 1C: 

Acquisitions: 

Restriction: 

OPTION 1Y: 

Acquisitions: 

Restriction: 

OPTION 2: 

Acquisitions: 


Restrictions: 


OPTION 2S: 

Acquisitions: 


D. Number of label able soybean dots must be >15. 

(Alternate option: 1C.) 

E. Number of label able N dots must be <7 or >15. 

(Alternate option: none.) 


Same as for option 1. 

Two-category labeling only: C and N (combine Y labels 

with N). 

Same as for option 1. 

Two-category labeling only: Y and N (combine C labels 

with N). 

A. At least one window Cl acquisition in Iowa, Illinois, or 
Nebraska; at least two window Cl acquisitions In other 
states. 

B. One window C2 acquisition. 

C. One or two window C3 acquisitions desirable but not 
requi red. 

A. Two-category labeling only: C and N. At least 15 dots 

in each category must be labeled. 

B. In areas where winter wheat Is planted following corn 

harvest (in excess of 5 percent of the pseudo-county), 
do not attempt this option. (Alternate option: none.) 

A. At least one window Cl acquisition In Iowa, Illinois, or 
Nebraska; at least two window Cl acquisitions In other 
states. 

B. One subwindow C2S1 acquisition. Select earliest 
acquisition available In the subwindow. 
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Restrictions: 

OPTION 3: 

Acquisitions: 

Restriction: 

OPTION 3S: 

Acquisitions: 


C. One subwindow C2S2 acquisition. Select latest 
acquisition available In the subwindow. 

D. One or two window C3 acquisitions desirable but not 
required. 

A. Number of label able corn dots must be >15. 

B. Number of label able N dots must be >15. 

(Alternate option: none.) 

A. At least one window Y1 acquisition in Iowa, Illinois, or 
Nebraska; at least two window Y1 acquisitions In other 
states. 

B. One window Y2 acquisition. 

C. One or two window Y3 acquisitions desirable but not 
requl red. 

Two-category labeling only: Y and N; 15 to 20 dots of each 

category required. 


A. At least one window Cl acquisition in Iowa, Illinois, or 
Nebraska; at least two window Cl acquisitions In other 
states. 

B. One subwindow C2S1 acquisition and one subwindow C2S2 
acquisition as in option 2S„ (Subwindow Y2S1 
acquisition may replace subwindow C2S2 acquisition If 
gap between windows C2 and Y2 Is <5 days.) 

C. One subwindow Y2S2 acquisition. 

D. One or two window C3 acquisitions desirable but not 
required. 
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Restrictions: 


Spring grains and sunflowers will be labeled N. 


OPTION 4S: 

Acquisitions: 

Restrictions: 

OPTION 5: 

Acquisitions: 

Restrictions: 


A. 

B. If the requirement for 15 to 20 dots cannot be met for 
three categories, combine the dots Into two categories 
and proceed as In less complex options. 

A. At least two window Y1 acquisitions. 

B. One subwindow Y2S2 acquisition. 

C. One or two window Y3 acquisitions desirable but not 
required. 

Label as categories Y and N; at least 15 dots must be 

labeled in each category. 


A. At least two window Y1 acquisitions. 

B. One window Y2 acquisition. 

C. One or two window Y3 acquisitions desirable but not 
required. 

The minimum of 15 Y and 15 N dots must be met. 
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SUMMER CROP LOGIC FLOW CHARTS 

The decision logic for labeling summer crops using the various labeling 

options is shown in figures G-l through G-5. An explanation of the 

terminology used in these flow charts is provided below. 

1. When a single window acquisition is referenced, the acquisition nearest 
the center of the window will be looked at first; the later acquisition 
will be meant in case of a tie. (For exceptions, see note 3 below.) 

2. When first, second, and third acquisitions are referenced, these should be 
taken to mean the acquisition nearest the center of the window, a later 
acquisition, and an earlier acquisition, respectively. 

3. In options 2S and 3S, reference to a subwindow C2S1 acquisition will be 
taken to mean the earliest acquisition available in the subwindow; 
reference to a subwindow C2S2 acquisition will be taken to mean the latest 
available acquisition, all other considerations being equal. The green 
number threshold of 13 on the C2S2 acquisition is tentative; further 
research is required to establish a relationship to the appropriate spring 
small grain growth stage and a corresponding corn growth stage. 

4. Option 3S includes option 2S in its entirety and requires manual input of 
Z dots (minus all corn dots) and N dots resulting from the option 2S 
spectral aids evaluation. Software to support processing of option 3S is 
not yet available; this option will be available at a later date. 
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Figure G-l.- Labeling logic for op 
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Figure G-2.- Labeling logic for option 2. An explanation of 
the terminology used in this flow chart is provided in 
notes 1 and 2 on page G-l. 
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Figure G-3.- Labeling logic for option 2S. An explanation of 
the terminology used In this flow chart Is provided In 
note 1, 2, and 3 on page G-l. 
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Figure G-4.- Labeling logic for option 3S . An explanation of 
the terminology used in this flow chart is provided in 
notes 1 through 4 on page G-l. 
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Figu r e G-5.- Labeling logic for option 4S . An explanation of 
the terminology used in this f’ow chart is provided in 
notes 1 and 2 on page G-l. 
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TECHNIQUE FOR USING PEAK GREENNESS TO DETERMINE 
PHENOLOGICAL MODEL CORRECTIONS 

! 

Prior methods available for determining the dates of labeling window occurrence 
(in order to select labeling acquisitions) have the following advantages and 
disadvantages: 


I. Analyst interpretation 


Advantage : Allows the analyst processing the segments to exert the 

most direct Influence on the selection of acquisitions. 


Disadvantages : • Requires several through-the-season acquisitions to 

allow the analyst to track the growth stages of each 
crop of interest and of other crops and confusion. 


• Depends on the analyst's prior training to identify 
the crop; poor acquisition selection will result from 
the analysts with lesser training. 

• Very time consuming. 


• Usually requires historical or current-year crop 

calendars (as an aid) for the geographic region being 
analyzed. 


2. Historical averages of growth stage occurrence 


Advantage : Allows complete automation of Window limit determination. 

Disadvantages : • Requires historical crop calendar data base. 

• Will not perform well in years having large deviations 
from the normal in growth stage occurrence, or in 
segments having large deviations. 

3. Thermal and photothermal phenol ogi cal models 

Advantages : • Allows growth stage predictions (and windows) to be 

completely automated. 
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• Adequately handles segment-to-segment variations for 
corn and soybeans. 


Pi sadvantages : • 


Depends on accurate temperature data from a 


sufficiently large network of reporting stations. 

Needs an accurate phenol ogical model. 

Does not adequately handle early and late seasons 
caused by rainfall fluctuations. 



The method of using peak greenness to determine labeling window midpoints, 
which are then used to correct phenol ogical models, has the following 
advantages and disadvantages: 


Advantages : • Allows results of skilled Interpretation of a few seg 

ments to be applied to many segments in an automated 
mode. 


• Same accuracy attainable with numerous acquisitions is 
attained in segments with minimal acquisitions. 

• Adequately handles segment-to-segment and season-to- 
season variations. 

• Complete crop calendar data base not essential. 

Disadvantages : • Depends on accurate temperature data from a network of 

reporting stations. 

• Interpretation of the crops and peak greenness 
determination of the few segments used to correct the 
models must be very accurate, since any errors made 
may be proliferated to other segments. 

• No segments can be processed until after the peak 
greenness stage of crops has been reached. 

• Either the 50-percent-planted date or the length of 
the growing season must be known. 
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Segments selected for using this technique must have several cloud-free 
acquisitions through the crop-growing season. This requirement Is necessary 
in order that curves of greenness through time can be constructed accurately 
enough to determine the peak greenness date for each selected crop pixel (dot) 
and to permit accurate temporal crop Interpretation. A planting acquisition 
and at least one acquisition before and after peak greenness are essential. 
Sufficient acquisitions must be present to allow construction of a curve for 
the earliest and latest fields. See figures H-l and H-2 for example curves of 
early and late developing crop dots. 

This technique also requires prior training in imagery Interpretation of the 
crop of interest. This training is essential, since the accuracy attained in 
window determination for a few segments will. In effect, be extended to many 
segments, because windows using this technique are used to adjust and train 
the models. 

After generating greenness through time plots for the segment, do the 
following: 

1. Using all of the cloud-free acquisitions for the segment, any available 
historical crop calendars, ground reports of conditions, and any other 
crop information for the segment year, make an accurate interpretation of 
a representative selection of as many pure corn and soybean dots as 
possible (15 is the desired minimum). 

2. Construct greenness through time curves of each dot and record the 
approximate date of occurrence of maximum greenness for each dot. This 
may be interpolated as a fraction of distance between two plotted dates. 

3. Calculate the mean of the peak greenness dates for each crop. 

4. Determine window limits as follows: 

• If both the 50-percent-planted date and the length of season are known: 

Cl = P - 20 through P„ + 21 
o o 

C2 = T - 15 through T + 15 
P P 
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Figure H-l.- Curve of early dot. 
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C3 = T + 0.36U ♦ 25 through T + 0.36L + 55 
p s P s 

Yl = P_ - 25 through P„ + 16 
o 3 o 

Y2 = T - 11 through T + 13 
P P 

Y3 = T + 0.37L + 10 through T + 0.37L e + 45 

P s p s 

where 

P 0 = date of 50 percent planted 
Tp = date of peak greenness 

L $ = length of season (I.e., number of days from 50 percent planted to 
50 percent harvested) 

NOTE: Neither window Cl nor window Yl may open earlier than day 110. 


• The alternate method to use when the length of season is known but the 
50-percent-pl anted date is not known is as follows: 

Cl = T p - 0.41L S - 25 through T p - 0.41l_ s + 16 

C2 = T - 15 through T + 15 

P P 

C3 = Tp + 0 . 3 6L s + 25 through T p + 0.36l_ s + 55 

Yl = T - 0.5L - 35 through T - 0.5L + 6 

p s P s 

Y2 = T - 11 through T + 13 
P P 

Y3 = T + 0.37L + 10 through T + 0.37L + 45 

p s P s 

where definitions of T p and L $ remain the same. 

NOTE: Neither window Cl nor window Yl may open earlier than day 110. 
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t The alternative method to use when the planting date, but not the 
length of season, Is known, Is as follows: 


Cl * P Q - 20 through P Q + 21 

C2 * T - 15 through T + 15 
P P 


C3 


■ °- 77 (- P TO—) + P 0 + 


30 through 0 




♦ + 60 
o 


Y1 • P - 25 through + 16 
o o 


Y2 * T p - 11 through T p + 13 


/T - P ft - 11\ 

V3 . 0.87(-i_| j 


+ P Q + 20 through 0 




+ P + 55 
o 


where the definitions of P Q and T p remain the same. 

NOTE: Neither window Cl nor window Y1 may open earlier than day 110. 


The second operation required in the technique Is to use the manually 
determined windows to correct the windows determined with the phenologlcal 
models. The necessary steps are as follows: 

1. Run the Cross-Zuber model and the Hajors-Johnson model to determine the 
model wIpjow midpoints: 

Cl = 50% planted 

C2 = 50% silked 

C3 = 50% mature 

Yl = 50% planted 

Y2 * 50% beginning pod fill 

Y3 = 50% mature 

2. Determine the midpoints of the windows determined using peak greenness. 
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3. Determine the window midpoint bias of each segment. Subtract the midpoint 
dates determined using peak greenness from the midpoint dates determined 
using the phenologlcal models. 

4. Determine the mean bias of each modeled window using the results of step 3. 

5. Change the sign of each window bias determined in step 4 and apply this 
number (of days) as a correction to each window output of the models. 
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APPENDIX I 

OUTPUT OF PROGRAM AIREPORT 
(SESSION 1) 

The semlautomated procedure presented In this document utilizes the Majors- 
Johnson and Cross-Zuber growth stage prediction models to establish the open- 
ing and closing dates of the time periods In which acquisitions are necessary 
for labeling corn and soybean crops. A sample printout of the Crop Report 
generated by this model Is shown on page 1-2. 

The Corn and Soybean Window Selection P.oort Is used by the analyst In the 
imagery screening session. This report contains the Initial acquisition 
selections within each window made by the computer. Ir addition, the options 
available for use at a particular latitude of a segment are printed in the 
report. These windows and subwindows, which are based on ^.he time periods 
predicted by the Majors-Johnson and Cross-Zubers growth stage prediction 
models, have been seasonally adjusted. 
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APPENDIX J 


PURE-PIXEL SELECTION USING MULTI TEMPORAL ACQUISITIONS 
AND EDGE GRADIENT ENHANCEMENT 

The method of selecting pure pixels is based on changes in the recorded Landsat 
digital values In channels 2 and 4. 

Each line of data is scanned, and the following values are recorded for each 
pixel : 

1. The difference between the value in a channel for each pixel and the 
following pixel 

2. The difference between the value in a channel for the same pixel and the 
pixel on the line below 

After recording the above differences for each pixel in both channels, the 
differences in each channel are histogrammed and differences below a 
predetermined threshold value are determined to be in the same field. 
Conversely, if a change is above the threshold in a channel, the two pixels 
being compared are deemed to be in different fields. The final decision is 
made as follows: If a pixel is below the threshold requirement in both 

channels, it is determined to be in the same field; if above the threshold in 
either channel or there is disagreement, the two pixels are determined to be in 
different fields. 

For multiple acquisitions, the above procedure is repeated for each acquisi- 
tion. If there is disagreement on whether pixels are in the same field, either 
in two-channel mode or across acquisitions, a pixel under examination is 
determined to be impure (in a different field). 
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APPENDIX K 


A COMPARISON OF BRIGHTNESS PROFILES AND 
GREENNESS-TO-BRIGHTNESS RATIO PROFILES 
FOR SEVERAL GRAINS AND GRASSES 

The graphic in this appendix illustrates how brightness, when combined with 
greenness as a ratio, can be a valuable aid in distinguishing between crops 
(such as wheat or oats) and pasture or grass hay, which can be confused with 
crops. 
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APPENDIX L 


EXAMPLES OF ANALYSIS AIDS 

This appendix contains examples of the spectral aids and listings used In the 
procedure session. The programs that produce these aids are specified In 
section 4.2 of this report. The examples contained In this appendix are as 
fol lows: 

1. Table A, page L-2: Green number-brightness listing for all acquisitions; 

ordered by dot number (418 dots) 

2. Table B, page L-3: Green number-brightness listing for all acquisitions; 

ordered by green number (418 dots) 

3. Table C, page L-4: Green number-brightness listing for the labeling 

acquisitions; ordered by labeling sequence (pure dots only) 

4. Unlabeled scatter plot, page L-5: Plot of all 418 dots 

5. Labeled scatter plot, page L-6; Plot of all pure summer crop (Z) dots 

6. Time trajectory plots, page L-7: Green number through time and brightness 

through time trajectories 

For products having multiple pages, only one page has been reproduced as an 
example. Example., of the total products are printed in reference 11. 
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APPENDIX M 


COMPUTER PROCESSING FUNCTIONS 

Section 4 of this report gives an end-to-end explanation of the processing 
functional flow of the charts contained In this appendix and how they relate 
to labeling and classification. Also, see reference 11. 
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Figure M-2.- Parameter estimation and dot labeling set creation. 
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Figure M-3.- Analysis aid production and final dot labeling. 
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Figure M-4.- Parameter estimation for classification 
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Figure M-5.- Classification and proportion estimation. 
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APPENDIX N 
TEST RESULTS 

OVERALL CS-4 VERIFICATION TEST RESULTS 
[SEVEN TEST SEGMENTS] 


Crop 

Classified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

37.15 

37.38 

-0.23 

-0.61 

Soybeans 

19.65 

20. 56 

-0.91 

-4.43 

Other 

43.20 

42.05 

♦1.15 

♦2.73 


DOT LABELING ACCURACY, X 


Ground 

truth 


Label 

C Y N~| 
C 86 7 7 


90 


1 


98 


Overal 1 

accuracy: 9 1.04 

Procedure error: 6.12 

Analyst error: 1.42 

Questionable 

ground truth: 1.42 


CLASSIFICATION VS LABELS, X 


CLASSIFICAT I ON VS GROUND TRUTH, % 


Classl flcatlon 


Classl flcatlon 


A 


Label 



C 

Y 

N 



C 

Y 

N 

C 

89 

8 

3 


C 

87 

8 

5 

Y 

9 

91 

0 

Dot 

ground 

truth 

Y 

9 

90 

1 

N 

8 

3 

89 


N 

5 

2 

93 


Purity of dots Initially selected by PURPIX program, X: 79.5 

Purity of labels after analyst acceptance, X: 95.8 
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CS-4 VERIFICATION TEST RESULTS 


Segment no. : 0144 
Year: 1979 
State: Iowa 


Crop 

Classified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

26.82 

19.88 

+6.94 

+34.91 

Soybeans 

16.41 

21.33 

-4.92 

-23.07 

Other 

56.76 

58.79 

-2.03 

-3.45 


Percentage of 
segment 

Inventoried: 96.8 


Number of dots 
(labels) In 
each category 
used to train 
classifier: 25 


i 


DOT LABELING ACCURACY. X 


Label 


C Y N 
C 100 0 0 


Ground 

truth 


93 0 


NO 0 100 


Overal 1 

accuracy: 97.3 

Procedure error: 1.35 

Analyst error: 0 

Questionable 

ground truth: 1.35 


Comment: Some hall dar.age and possibly poor stands caused some soybeans to be classified 

as corn. 
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Cl ass 1 flcatlon Classification 
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N ~ 
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C 

88 

4 

8 

Dot 

C 

87 

4 

9 

Label 

Y 

12 

88 

0 

qround 

truth 

Y 

12 

88 

0 


N 

4 

4 

92 


N 

4 

4 

92 


Purity of dots Initially selected by PllRPIX program, X: 88.9 

Purity of labels after analyst acceptance, X: 92.0 
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CS-4 VERIFICATION TEST RESULTS 


Segment no.: 0145 
Year: 1979 
State: Iowa 


Percentage of 
segment 

Inventoried: 96.1 


Number of dots 
(labels) In 
each category 
used to train 
classifier: 15 


DOT LABELING ACCURACY. X 


Crop 

Classified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

29.51 

32.36 

-2.85 

-8.81 

Soybeans 

24.11 

16.63 

♦7.48 

+45.0 

Other 

46.37 

51.01 

-4.64 

-9.00 


Label 

C Y N 
C 100 0 0 


Ground 

truth 


83 


0 


N 0 0 100 


Overal 1 

accuracy: 93.3 

Procedure error: 6.7 

Analyst error: 0 

Questionable 
ground truth: 0 


Comment: Small amounts of corn classified as soybeans had a high relative effect on the 

total classified soybeans. 
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Classl flcatlon Classification 



C Y N 


C Y N 


C 73 27 0 

Dot 

C 92 8 0 

Label 

Y 13 87 0 

ground 

truth 

Y 11 89 0 


_N 6.5 6.5 87 


N 6.5 65 87 


Purity of dots Initially selected by PURPIX prc^ram, %: 73. 5 

Purity of labels after analyst acceptance, X: 100 
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CS-4 VERIFICATION TEST RESULTS 


Segment no.: 0856 
Year: 1979 

State: Indiana 


Crop 

Clarified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

27.28 

31.42 

-4.14 

-13.22 

Soybeans 

28.00 

25.77 

♦2.23 

-8.70 

Other 

44.72 

42.81 

♦1.91 

♦4.50 


Percentage of 
segment 

Inventoried: 96.9 


Number of dot* 

(labels) In 
each category 
used to train 
classifier: 25 


DOT LABELING ACCURACY. X 


Label 

C Y n' 

C 85 11 4 


Ground 

truth 


8 92 


0 


N 


0 


0 100 


Overal 1 

accuracy: 92.0 

Procedure error: 4.0 

Analyst error: 1.3 

Questionable 

ground truth: 2.7 


Comment: Some late corn fields had spectral characteristics of soybeans. 


CLASSIFICATION VS LABELS. X 


CLASSIFICATION VS GROUND TRUTH. % 


Classification Classification 
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Dot 
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0 

Y 

0 

100 

0 

ground 

trutTT 
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8 

92 

0 

N 

4 

0 

96 
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0 

0 

100 


Purity of dots Initially selected by PURPIX program, X: 78.9 

Purity of labels after analyst acceptance, X: 98.7 
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Segment no. : 0864 
Year: 1978 
State: Iowa 


Crop 

Classified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

42.47 

46.21 

-3.74 

-8.09 

Soybeans 

10.61 

11.95 

-1.34 

-11.21 

Other 

46.92 

41.84 

+5.08 

♦12.21 


Percentage of 
segment 

Inventoried: 97.) 


Number of dots 
(labels) In 
each category 
used to train 
classifier: 17 


DOT LABELING ACCURACY. X 


Labe l 
C Y 


C 89 


N 

11 


Ground 

truth 


94 


0 10G 


Overal 1 

accuracy: 94.1 

Procedure error: 3.9 

Analyst error: 2.0 

Questionable 
ground truth: 0 


Comment: Some early maturing corn *ields were below corn thresholds and were classified 
as "other"; some soybean fields were called corn (possibly poor stands). 
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Classification 
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C 

32 

12 
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83 

6 

11 






Dot 





Ubel_ 

Y 

6 

94 

0 

ground 

Truiir 

Y 

6 

94 

0 


N 

6 

0 

94 


N 

0 

0 

100_ 


Purity of dots initially selected by C URPIX program, X: 73.2 

Purity of labels after ar alyst acceptance, X: 96.1 
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CS-4 VERIFICATION TEST RESULTS 


Segment no.: 0881 
Year: 1978 
State: Io*a 


Crop 

Classified, 

X 

^ound 
truth, X 

Error, 

X 

Relative 
error, X 

Corn 

42.11 

44.33 

-2.22 

-5.00 

Soybeans 

10.53 

8.00 

+2.53 

+31.63 

Other 

47.36 

47.67 

-0.31 

-0.65 


Percentage of 
segment 

Inventoried: 98.04 


Number of dots 
(labels) In 
each category 
used to train 
classifier: 18 


DOT LABELING ACCURACY. X 


Label 

C Y N 
C 75 21 4 


Ground 

truth 


0 100 


0 


N 0 6 94 


Overal 1 

accuracy: 86.8 

Procedure error: 9.4 

Analyst error: 3.8 

Questionable 

ground truth: 0 


Comment: Some late corn fields still had spectral characteristics of soybeans. 
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Classification 


Cl assi flcation 
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94 

6 
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C 

79 

17 
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Dot 





Y 

17 

83 

0 

ground 

truth 

Y 

0 

100 

0 

N 

11 

0 

89 


N 

12 

0 

88_ 


Purity of dots Initially selected by PURPIX program, X: 71.2 

Purity of labels after analyst acceptance, X: 92.6 
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Segment no. : 0877 
Year: 1978 

State: Iowa 


Percentage of 
segment 

Inventoried: 79.3 


Number of dots 
(labels) In 
each category 
used to train 
classifier: 25 


DOT LABELING ACCURACY. X 


Crop 

Classified, 

% 

Ground 
truth. % 

Error, 

t 

Relative 
error, l 

Corn 

46.95 

48.30 

-1.35 

-2.80 

Soybeans 

20.90 

24.88 

-3.98 

-15.68 

Other 

32.15 

26.82 

♦5.33 

+19.87 


Label 


C 


Y 



C 78 6 16 


Ground 

truth 


0 90 10 


NO 0 100 


Overal 1 

accuracy: 86.8 

Procedure error: 8.8 

Analyst error: 1.5 

Questionable 

ground truth: 2.9 


Comment: Variety of factors (Including early maturity) caused small amount of corn to be 

called "other"; some soybeans labeled "other" — failure to detect planting. 


CLASSIFICATION VS LABELS, % 
Classification 


Label 


CLASSIFICATION VS GROUND TRUTH, % 
Classi f 1cation 
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C 

100 

(. 

0 


c 

91 

3 

6 
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12 

88 

0 

Dot 

ground 

truth 
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0 

95 

5 

N 



12 

4 

84 

— J 


N 

0 

0 

100 


Purity of dots initially selected by PURP!X program, %: 82._2 

Purity of labels after analyst acceptance, l: 96.0 
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CS-4 VERIFICATION TEST RESULTS 


Segment no. : 0883 
Year: 1979 
State: Iowa 


Crop 

Classified, 

X 

Ground 
truth, X 

Error, 

X 

Relative 
error, % 

Corn 

44.92 

39.17 

♦5.75 

+14.7 

Soybeans 

26.98 

35.39 

-8.41 

-23.76 

Other 

i 

28.10 

25.44 

+2.66 

+10.5 


Percentage of 

segment 

Inventoried: 

Number of dots 
(labels) In 
each category 
used to train 
classifier: 


95.4 


•k 

19 


DOT LABELING ACCURACY. X 


Label 

C Y N ' 
C 88 0 12 


Ground 

truth 


78 


0 


N 


0 


6 94 


Overal 1 

accuracy: 86.0 

Procedure error: 10.5 

Analyst error: 3.5 

Questionable 

ground truth: 0 


Comment: Extensive hall damage caused corn to be called "other" and soybeans to be called 

corn. 


CLASSIFICATION VS LABELS. % 


CLASSIFICATION VS GROUND TRUTH, % 


Label 


Classl flcatlon Classification 
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Y 

N 



C 

Y 

N 

C 

84 

11 

5 


C 

94 

0 

6 





Dot 





Y 

r 

95 

0 

ground 

truITT 

Y 

13 

87 

0 

N 

21 

0 

79 _ 


N 

l7 

0 

83 


Purity of dots Initially selected by PURPIX program’ %: 85. 7 

Purity of labels after analyst acceptance, X: 96.5 
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